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COMPUTED  AND  EXPERi:4ENTAL  HUGONIOTS  FOR  UNREACTED  POROUS  HIGH 
EXPLOSIVES 

The  work  described  in  this  report  was  carried  out  under  Task 
MAT-03L-000/R011-01-01 ,  IR  159,  Transition  from  Deflagration  to 
Detonation  of  NAVSURWPNCEN  Independent  Research  Program.  Both 
computed  and  experimentally  measured  dynamic  compression  curves 
(Hugoniot  curves)  for  unreacted  TNT  are  presented.  Also,  our 
computations  for  several  other  porous,  unreacted  explosives  are 
compared  with  experimental  data  from  the  Russian  literature.  The 
work  shows  that  Hugoniots  can  be  computed  for  unreacted,  porous 
explosives.  This  result  is  important  because  it  simplifies  the 
determination  of  the  sensitivity  of  porous,  granular  explosives  to 
shock  stimuli. 
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INTRODUCTION 


A  knowledge  of  the  Hugoniot  of  unreacted  explosives  is  required 
to  interpret  transitional  phenomena  in  explosives  and  to  determine 
the  shock  sensitivity  of  explosives^.  For  example,  the  Naval  Surface 
Weapons  Center  Large  Scale  Gap  Test  (LSGT) ^ ^  is  used  to  determine 
the  relative  shock  sensitivity  of  explosives.  The  LSGT  is  calibrated 
to  give  the  pressure  in  the  gap  material  as  a  function  of  gap  length. 
But  after  a  determination  of  the  critical  gap  pressure  for  initiating 
detonation,  the  question  still  remains:  What  is  the  critical  pressure 
transmitted  to  the  test  explosive?  This  question  cannot  be  answered 
without  a  knowledge  of  the  Hugoniot  of  the  unreacted  explosive. 

The  usual  method  employed  to  obtain  the  Hugoniot  of  a  material  is 
to  perform  a  series  of  experiments  which  measure  two  of  the  following 
variables:  shock  velocity,  particle  velocity,  pressure,  and  volume. 
There  are  two  drawbacks  to  this  approach  for  porous  (granular) 
explosives.  (1)  Because  porous  explosives  react  very  readily,  it  is 
difficult  to  measure  any  of  the  above  variables  for  the  unreacted 
state.  (2)  Because  a  new  Hugoniot  is  needed  for  every  density,  this 
approach  rapidly  becomes  expensive.  However,  computational  methods 
have  been  developed  which  yield  the  Hugoniot  of  inert,  porous 
materials  when  the  Hugoniot  of  the  non-porous  material  is  known.  If 
these  computational  methods  can  applied  to  porous  explosives,  much 
work  and  expense  can  be  saved. 

In  this  study,  the  authors  have  applied  these  computational 
methods  to  porous  TNT  for  a  ran-^e  of  densities  from  1.60  to  0.98  g/cc. 
This  method  was  confirmed  by  experimentally  obtaining  the  Hugoniot 
for  unreacted  TNT  at  a  density  of  0.98  g/cc  over  the  range  of  particle 
velocities  of  0.350  -  1.720  mm/ysec  and  of  pressures  of  3. 6 ^  55 . 0  kbar . 
The  electromagnetic  velocity  (EMV)  gage**  was  used  in  conjunction  with 
the  known  Hugoniot  data  of  the  attenuator  material  (PMMA)  to  obtain 
the  Hugoniot  data  for  unreacted  TNT  at  0.98  g/cc.  The  EMV  gage 
follows  particle  velocity  as  a  function  of  time.  The  gage  was  mounted 
at  the  inert/TNT  interface  and  responded  to  the  particle  velocity  at 
the  Interface.  In  the  pressure  range  (3.6  -  55.0  kbar)  covered  in 
this  work,  some  finite  induction  period  is  required  before  reaction 
occurs  in  the  TNT.  The  EMV  gage  leads  were  covered  by  a  thin  layer 
of  mica  to  insulate  them  from  any  reaction  products. 
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COMPUTATIONAL  MODELS 

A  computational  model  which  reproduces  the  response  of  porous 
materials  to  shock  has  been  used  by  several  researchers^'®'^.  It  has 
been  usefully  applied  to  a  wide  variety  of  materials  —  metals, 
graphite,  plastics,  boron  nitride  and  others  —  all  of  which  are  non¬ 
reactive.  Successful  application  of  this  model  to  potentially  reactive 
materials  such  as  porous  explosives  has  the  advantage  of  expediting 
the  studies  of  these  materials  by  saving  much  experimental  effort. 

The  model  requires  a  knowledge  of  the  equation  of  state  of  the  material 
under  study  in  solid  (voidless)  form.  This  is  furnished  by  a  curve  or 
relation  giving  the  Hugoniot  of  the  solid  and  a  relation  for  computing 
the  Gruneisen  parameter,  r (v) .  Because  of  the  lack  of  data,  we  have 
had  to  assume  that  the  Gruneisen  parameter  is  a  constant  in  our  com¬ 
putations. 

The  relation  between  the  solid  Hugoniot  and  the  porous  Hugoniot 
is  derived  as  follows®.  The  Gruneisen  equation  of  state  is  written  as 

Pf  -  Ps  =  (r/v)  (Ef  -  Eg)  (1) 

where  P  is  the  pressure,  v  is  the  specific  volume,  E  is  the  specific 
internal  energy,  and  the  subscripts  f  and  s  refer  to  quantities  in  the 
porous  and  solid  material  respectively  at  the  same  volume  v.  Because 
the  shear  strengths  of  explosives  are  not  known,  we  use  the  pressure 
to  approximate  the  principal  stress  in  the  direction  of  shock  propaga¬ 
tion.  The  Rankine-Hugoniot  energy  equations  for  both  materials  are 


Ef  -  Eof  =  Pf(Vof  -  v)/2 


-  Eos  =  Ps<Vos 


-  v)/2 


where  the  added  subscript,  o,  refers  to  the  initial,  or  zero  pressure, 
state  of  the  material.  The  assumption  is  made  that  Eof  =  ^og;  Boade® 
and  others®  justify  this  assumption  on  the  basis  that  the  two  energies 
differ  only  bv  the  additional  small  amount  of  surface  energy  in  the  porous 
material.  With  this  assumption  and  substitution  of  equations  (2)  and 
(3)  into  equation  (1) ,  we  obtain 

(vos  -  V  -  2v/r| 

P^:  =  P.-  _  (4) 


Vof  -  V  -  2v/r  j 


For  some  materials,  the  relation  between  the  shock  velocity,  U, 
and  the  particle  velocity,  u,  is  linear,  i.e., 

U  =  Co  +  S  u,  (5) 

where  Cq  is  the  sound  speed  at  zero  pressure  and  S  is  the  slope  of  the 
U,u  curve.  For  such  materials,  a  convenient  representation  for  Pg  can 
be  derived  by  using  the  Rankine-Hugoniot  equations 

P  -  =  u  U  /  and 

U  -  Uq  =  (U  -  u)  Vq/  V.  (7) 
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Solving  for  the  pressure,  there  results 


(Vqs  -  V) 

[Vos  -  S(Vos  -  v)]2 


(8) 


pQ  (  1  bar)  is  considered  negligible  compared  to  P  in  Eq.  (6)  and  Uq 
in  Eq.  (7)  is  zero  in  applications  considered  here.  With  these 
approximations,  equations  (4)  and  (8)  will  yield  the  Hugoniot  for  the 
porous  material. 


Equations  (4)  and  (8)  are  included  in  a  computer  code  written  by 
Herrmann^.  (Also  included  is  a  method  for  calculating  the  behavior  of 
porous  materials  while  the  pores  are  collapsing  —  the  crushup  process. 
This  method  is  generally  known  as  the  P-a  model  where  a  is  a  parameter 
which  describes  the  compaction  process.  The  required  data  for  this 
code  are:  the  Hugoniot  of  the  solid  material,  a  value  for  r (v)  (we 
assume  that  P  is  a  constant  in  this  report)  ,  the  specific  vol\ime  of 
the  porous  material,  the  crushup  pressure  (Pq) ,  and  the  pressure  at 
the  elastic  limit  for  the  porous  material) .  A  computer  code  is 
described  in  Appendix  A  which  solves  Eqs.  (4)  and  (8),  i.e..  it  does 
not  treat  the  compaction  process. 

Throughout  this  report  we  will  be  exploiting  only  a  small  portion 
of  the  work  published  by  Herrmann’.  As  noted  above,  we  do  not  use  his 
scheme  for  calculating  the  pressure  in  the  crushup  region.  Neither 
do  we  use  his  more  complicated  relations  for  describing  the  Hugoniot 
of  the  solid  material.  These  latter  relations  remove  the  restrictions 
of  linearity  of  the  relationship  between  shock  and  particle  velocities. 


Although  the  above  treatment  seems  to  be  the  currently  most 
successful  treatment  of  granular  materials,  other  methods  of  determin¬ 
ing  the  Hugoniot  are  also  used.  In  particular,  the  Russian  literature 
on  explosives  frequently  cites  Voskobinikov  et  al®.  Starting  with  a 
generalized  Hugoniot  for  organic  liquids,  these  workers  derived  an 
expression  for  ar.  organic  liquid  containing  bubbles.  They  considered 
it  applicable  to  porous  solids  as  well.  Their  final  expression  is 


/ 

= 


1  + 


2P 


(Y+1)  u2 


(vjf  -  Vq3) 


(9) 


where  y  =  (Cp/Cy)  for  air  and  Vof  is  given  by 


-St  =  ''of  +  ^  ^  ''oair 


(10) 


where  <5  is  the  mass  fraction  of  the  air  in  the  porous  material. 

In  a  recently  translated  article,  Dremin,  Shvedov,  and  Avdonin^ ° 
use  a  more  involved  method  for  calculating  the  Hugoniots  of  porous 
materials.  The  method  is  described  in  some  detail  by  Al * tshuler  et  al 
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in  a  paper^*  published  in  Russia  in  1958.  Zel*dovich  and  Raizer'^also 
describe  the  method.  The  method  starts  with  the  use  of  the  caloric 
equation  of  state 

^  ^  ST  =  -  Ee)/v  (11) 

where  the  subscript  c  refers  to  the  0'’K  isotherm.  This  is  the  chief 
difference  between  this  method  and  the  one  used  by  Herrmann  and  others 
who  used  the  Hugoniot  of  the  solid  as  the  reference  thermodynaraic  path. 
Equation  (11)  is  specialized  to  describe  a  Hugoniot  by  replacing  E  by 
0.5  P  (vq-v)  ,  the  change  in  energy  across  a  shock.  There  result-  ; 

V  +  TEp  =  (Vq  -  V  -  2v/r)  .  (12) 

dv  ^2 

The  subscript  H  means  that  the  pressure  is  evaluated  along  the  Hugoniot 
curve.  Al*tshuler  et  al  show  how  Eq.  (12^  can  be  converted  into  an 
integral  equation  which  can  be  solved  nume/icallv .  The  reason  for 
doing  this  is  that  is  a  function  of  the  volume,  which  has,  so  far, 
not  been  specified.  In  order  to  solve  the  equation,  the  function 
?h(v)  must  be  known,  usually  from  experiments.  If  the  variable  T  is 
also  treated  as  a  function  of  the  volume?  it  also  must  be  specified. 
Dremin  et  al  treated  F  as  a  constant,  as  will  be  done  throughout  this 
discussion. 

Once  Ec (v)  is  known,  the  Hugoniot  of  the  same  material  in  porous 
form  can  be  obtained  from  Eq,  (12)  .  That  is, 

2v 

PH(''of.v)  =  (  ^  -^  +  2Ec)/(v^f  -  V  -  2v/r),  (13) 

where  v  £  is  the  initial  specific  volume  of  the  porous  material.  This 
comes  about  because  in  the  derivation  of  Eq,  (12),  no  assumption  was 
made  as  to  whether  the  material  was  continuous  or  discontinuous.  As 
in  the  derivation  of  Eq.  (4),  it  is  assumed  here  that  the  surface 
energy  of  the  porous  material  can  be  neglected. 

Because  Ec(v)  is  known  as  a  function  of  v,  Eq.  (13)  can  also  be 
solved  numerically,  giving  ^  function  of  v  for  any 

specified  value  of  v  £.  This  Dremin,  Shvedov,  and  Avdonin have  done 
for  RDX  (1.0  g/cc) ,  tetryl  (0.B6  q/cc]  ,  PETN  (0.82  g/cc)  and  ammonium 
nitrate  (0.86  g/cc) .  Their  results  agree  well  with  their  experimental 
data  as  is  shown  in  a  later  section  of  this  report.  V7e  used  the 
simpler  model,  Eqs .  (4)  and  (8),  to  compute  Hugoniots  for  the  same 
explosives.  These  results  will  be  discussed  in  a  later  section. 

In  concluding  this  section  of  the  report,  we  must  point  out  that 
the  above  treatments  of  shock  compression  of  porous  materials  neglects 
what  is  Called  the  "crushup”  process.  This  is  the  process  which 
closes  the  voids  in  the  porous  media.  Herrmann  has  extended  his  work 
to  compute  flow  variables  during  crushup.  We  do  not  use  that  part  of 
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his  work  because  we  lack  the  necessary  input  data.  Required  data  are 
the  sound  speed  in  the  porous  material,  its  elastic  limit,  and  the 
pressure  at  which  it  becomes  void  free.  More  elaborate  models  of  the 
crushup  process  are  available,  which,  because  they  require  even  more 
data,  cannot  at  present,  be  applied  to  our  work.  As  for  the  work  of 
Dremin  et  al,  they  warn  about  extrapolating  their  data  to  pressures 
below  10  kbar  because  of  the  possibility  of  splitting  of  the  shock 
waves.  Hence,  they  are  aware  of  the  phenomena  of  crushup.  Based  on 
the  translated  articles  we  have  seen,  they  do  not  have  a  computational 
model  which  accounts  for  the  crushup  process.  In  one  Russian  article, 
it  is  stated:  "There  is  still  no  reliable  description  of  the  behavior 
of  a  porous  material  dynamically  compressed  to  a  density  less  than  the 
density  of  the  single  crystal."^® 

In  the  next  section  of  this  report,  the  experimental  method  used 
to  obtain  Hugoniot  data  for  0.98  g/cc  TNT  is  described.  In  a  later 
section,  the  experimental  data  are  compared  with  results  computed 
using  seme  of  the  models  described  above. 


EXPERIMENTAL 

The  experimental  configuration  used  in  this  work  is  shown  in 
Figure  1.  The  pentolite  booster  and  polymethyl  methacrylate  (PMMA) 
cylinder  constitute  the  pressure-donor  and  pressure-attenuator  parts, 
respectively,  of  the  NAVSURWPNCEN  Large  Scale  Gap  Test  (LSGT)^'^. 

The  50/50  pentolite  booster  consists  of  two  pellets,  50.8  mm  diameter 
by  25.4  mm  long,  at  a  density  of  1.56  g/cc.  The  PMMA  cylinder  is 
50.8  mm  diameter  and  is  machined  to  the  desired  length,  x. 

The  electromagnetic  velocity  (EMV)  gage  consists  of  a  rectan¬ 
gular  loop  of  aluminum  foil  0.025  mm  thick  and  *-2  mm  wide;  the 
effective  gage  base  length^**  used  in  this  work  was  .1.11  cm.  The  EMV 
gage  used  in  this  work  is  in  the  radial  configuration  (i.  e.,  the 
leads  come  out  the  side  of  the  charge)  because  this  configuration 
facilitates  the  assembly  of  the  setup.  The  EMV  gage  is  mounted  on 
the  surface  of  the  PMMA  cylinder  with  chloroform  and  the  TNT  is  then 
hand  pressed  to  the  desired  density  of  0.98  ±  0.03  g/cc.  The  TNT 
(Lot  X517)  is  a  coarse,  granular  material  and  has  an  average  particle 
size  of  approximately  400y.  This  particle  size  was  estimated  by  a 
partial  sieve  analysis. 

The  experimental  assembly  is  placed  in  a  magnetic  field  of  "400 
gauss.  This  field  is  generated  by  a  pair  of  Helmholtz  coils  which 
have  a  35.5  cm  diameter  and  750  turns/coil  of  #14  wire.  The  coils 
are  energized  by  a  20  amp,  120  7olt  constant  current  supply.  The 
instrumentation  and  data  reduction  method  used  in  this  work  are 
described  in  Reference  (4) . 

The  pressure-donor  and  pressure-attenuator  part  of  the  LSGT  was 
chosen  for  this  work  because  the  LSGT  is  a  calibrated  pressure  vs 
distance  (P,x)  and  particle  velocity  vs  distance  (u,x)  system  which 
can  be  used  with  Helmholtz  coils  without  destroying  them.  The  main 
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FIG.  1  EXPERIMENTAL  ARRANGEMENT 
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drawbacks  of  this  system  are;  (1)  the  shock  front  in  the  PMMA  is 
curved,  and  (2)  PMMA  is  a  viscoelastic  solid.  Earlier  work  with  the 
EMV  gage**  in  curved  shock  fronts  showed  that  the  main  effect  of  the 
curved  fronts  was  to  lengthen  the  rise  time  of  the  output  signals. 


ANALYSIS  OF  DATA 

Electromagnetic  Velocity  Gage  Technique 

The  EMV  gage  record  yields  a  particle  velocity  vs  time  profile 
of  the  shock  loaded  TNT  at  the  PMMA/TNT  interface.  The  governing 
equation  for  this  method  is 

u(t)  =  V(t)/TH  10"")  (14) 

where  u(t)  is  the  particle  velocity  as  a  function  of  time  in  mm/psec, 
V',t)  is  the  output  voltage  as  a  function  of  time  in  volts,  H  is  the 
magnetic  field  in  gauss,  and  5,*  is  the  effective  gage  base  length'** 
in  mm.  For  the  radial  configuration,  Z*  is  the  distance  between  the 
middle  of  the  leads,  see  Figure  1. 

The  u,t  data  are  fitted,  for  most  cases,  to  a  straight  line  by 
least  squares  using  the  data  from  the  peak  recorded  value  of  u  to  a 
point  0.5  ysec  beyond  this  value.  This  is  the  same  approach  as  was 
used  in  Reference  3.  The  velocity  intercept,  Ug,  is  then  assumed  to 
be  the  particle  velocity  at  the  shock  front.  For  the  lowest  pressure, 
the  u,t  curve  had  a  long  risetime  ("'500  ns)  and  did  not  have  a  sharp 
peak.  For  this  record,  we  simply  use  the  peak  of  the  observed  curve 
for  the  shock  induced  particle  velocity. 

The  risetime  of  the  u,t  profiles  in  the  porous  TNT  is  longer  than 
that  obtained  in  a  voidless  solid  by  a  factor  of  two  or  more.  This  is 
apparently  due  to  the  time -dependent  collapse  of  the  pores'®  . 

Experiments  are  performed  using  attenuators  of  different  length 
so  that  we  have  several  pairs  of  values  of  Ug  and  x,  the  length  of 
the  PMMA  attenuator.  From  these  data,  and  the  known  Hugoniot  of  PMMA 
we  can  obtain  the  pressure,  P^,  at  the  interface.  The  process  for 
doing  this  is  discussed  in  the  section  "Impedance  Matching". 

Effects  of  Reaction  on  Gage  Output 

Reaction  is  always  a  potential  problem  in  the  experimental  deter¬ 
mination  of  the  unreacted  Hugoniot  of  a  shock  sensitive  material.  In 
this  study,  the  experimental  problem  created  by  reaction  is  that 
explosives,  and  especially  TNT,  become  electrically  conductive  and 
affect  the  gage  output.  The  TNT  is  going  to  react  to  some  degree, 
during  the  0-5ysec  observation  time,  especially  at  high  pressures. 
However,  even  if  the  TNT  starts  to  react,  the  initial  value  of  the 
particle  velocity  at  the  interface  should  lie  on  the  unreacted  Hugoniot. 
The  reason  for  this  is  the  following.  The  shock  emerging  from  the 
PMMA  is  non-reactive.  Any  reaction  which  occurs  at  the  interface. 
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where  the  gage  is  located,  must  start  at  zero  reaction  and  follow  some 
type  of  rate  law.  Thus  the  initial  shock  state,  and  the  particle 
velocity  associated  with  this  state,  must  be  on  the  unreacted  Hugoniot. 
The  problem  is  to  obtain  this  initial  value  without  interference  from 
reaction. 

Before  this  problem  can  be  solved,  some  idea  of  how  reaction 
affects  the  gage  is  necessary.  The  EMV  gage  has  been  used  in  a  non¬ 
reactive  material,  PMMA,  in  an  earlier  study^.  Figure  2  is  a  typical 
gage  record  when  the  gage  is  located  in  PMMA  at  a  distance,  x  =  25  mm, 
from  the  booster.  Notice  the  fairly  smooth  falloff  of  particle 
velocity  with  time.  Since  PMMA  is  the  inert  used  in  this  work  the 
particle  velocity  vs  time  profile  recorded  at  the  PMMA/TNT  interface 
should  resemble  Figure  2  if  the  TNT  reaction  does  not  affect  the  gage 
output.  In  order  to  see  how  reaction  affects  the  gage  output,  an  EMV 
gage  was  placed  in  the  TNT  at  a  position  where  the  TNT  was  reacting. 
Figure  3  shows  the  results.  The  electrical  conduct!’,  ity  associated 
with  the  reaction  products  of  the  TNT  apparently  sliorts  out  the  gage 
leads  and  causes  the  output  of  the  gage  to  decrease  rapidly.  The 
initial  particle  velocity,  which  is  the  velocity  at  time  t  =  0,  is 
the  quantity  of  interest.  To  obtain  this  velocity,  the  u,t  profiles 
must  be  extrapolated  back  to  t  =  0.  The  extrapolated  value  using  a 
u,t  profile  similar  to  Figure  2  is  more  reliable  than  the  value 
obtained  using  a  u,t  profile  similar  to  Figure  3.  If  the  gage  leads 
can  be  adequately  insulated,  the  u,t  profiles  obtained  from  the  EMV 
gage  at  the  PMMA/TNT  interface  should  resemble  the  u,t  profile  of 
Figure  2. 

To  circumvent  the  problem  of  conductivity,  especially  at  high 
pressures,  it  was  decided  to  cover  the  EMV  gage  leads  with  an  insulator. 
A  thin  (0.001  in.)  layer  of  mica  was  used.  A  shot  was  fired  with  two 
gages  at  the  interface,  (where  x  =  12.7  mm,  the  length  of  the  PMMA 
attenuator)  one  with  leads  covered  with  mica,  the  other  bare.  Figure  4 
shows  the  resulting  u,t  profiles.  The  lower  curve  is  the  bare  gage 
result.  It  is  very  similar  to  the  record  from  the  gage  located  in 
the  reacting  TNT,  see  Figure  3.  The  t  =  0  value  of  the  particle 
velocity  obtained  using  this  curve  would  not  be  very  reliable.  The 
upper  curve  is  the  mica  covered  gage  result.  This  curve  is  very 
similar  to  the  u,t  profile  in  PMMA,  see  Figure  2.  Thus,  covering  the 
gage  leads  with  an  insulator  (mica)  prevents  conductivity  from 
affecting  the  gage  output  up  to  at  least  55  kbar  which  was  the  highest 
pressure  used  in  this  study. 

Impedance  Matching 

When  a  shock  encounters  an  interface  between  materials  having 
different  shock  impedances,  a  shock  is  transmitted  across  the  inter¬ 
face,  and  a  wave  is  reflected  back  into  the  originally  shocked  material . 
Here  impedance  means  the  product  of  the  density  and  the  shock  velocity. 
The  reflected  wave  will  be  a  rarefaction  if  the  impedance  of  the 
second  material  is  less  than  that  of  the  first.  Otherwise,  it  will 
be  a  shock  wave.  Solving  for  the  strengths  of  the  reflected  and 
transmitted  wave  can  be  done  graphically,  or  analytically,  and  is 
called  impedance  matching  (or  mirror  image)  technique. 
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We  have  the  Hugoniot  curve,  i.e.,  a  relation  between  P  and  u,  for 
PMMA  for  forward  facing  shocks.  This  was  derived  from  the  relations^ 

U  =  2.561  +  1.595U  u  i  0.536, 

and 

U  =  2.7228  +  4.0667U  -  lO.SOSlu^  (15) 

+  10.691u3,  u  <  0.536, 


where  U  and  u  are  the  shock  velocity  and  particle  velocity,  respectively 
in  mm/ysec.  Substituting  equation  (15)  into  equation  (6)  with  Pq 
assumed  to  be  zero  leads  to 

P  =  10Pq(2.561u  +  1.595u2)  u  >  0.536, 

and 

P  =  10pj^(2.7228u  +  4.0667u2  -  10.9051u3  (16) 

+  10.691u‘*)  ,  u  <  0.536, 


where  P  is  in  kbars  and  Pq  is  the  density  of  the  PMMA,  which  is 
1.185  g/cc.  Equations  (16)  are  represented  by  the  curve  labeled  1 
in  Figure  5  which  represents  only  the  first  quadrant  of  the  P,u 
diagram.  This  curve  is  for  the  case  of  shocks  which  are  propagating 
in  the  forward,  or  +x,  direction  in  a  conventional  t,x  diagram.  Had 
the  shocks  been  backward  facing,  i.e.,  propagating  into  PMMA  in  the 
-X  direction,  the  curve  would  pass  through  the  origin  of  Figure  5 

with  negative  slope,  i.e.,  it  would  lie  in  the  second  quadrant  of  the 

P,u  diagram.  Such  a  curve  can  be  translated  to  the  right,  see  the 
curve  labeled  2  in  Figure  5.  This  is  the  P,u  curve  for  backward 
facing  shocks  in  PMMA  which  had  a  particle  velocity  of  2ui  in  the  for¬ 
ward  direction  prior  to  being  shocked.  Curve  2  is  said  to  be  the 
reflection  of  curve  1  about  the  line  u  =  u^.  Mathematically,  any 
such  curve  is  described  by  replacing  u  in  equation  (16)  by  (2u2-u)  . 

As  noted  above,  the  wave  reflected  back  into  the  PMMA  may  be  a 

shock  or  a  rarefaction.  If  it  is  a  shock,  values  of  P  and  u  induced 
by  it  will  lie  on  the  part  of  curve  2  above  the  intersection,  point  A 
in  Figure  5.  This  part  of  curve  2  is  exact  for  reflected  shocks.  If 
the  reflected  wave  is  a  rarefaction,  we  use  the  part  of  curve  2  which 
lies  below  point  A.  But  flow  variables  in  a  relief  wave  (simple  wave) 
are  governed  by  ad-'abatic  or  isentropic  relations.  Curve  2  is  exact 
only  for  shock  compression,  a  process  which  produces  entropy.  To 
solve  for  a  relief  wave  in  the  PMMA,  the  curves  should  be  those  for 
an  adiabatic  process. 


A  common  approximation  is  made  by  noting  that  the  transition 
through  a  simple  wave  gives  values  of  the  flow  variables  which  agree 
to  second  order  with  those  produced  by  a  weak  shock .  So  to  second 
order,  the  segment  of  curve  2  below  point  A  of  Figure  5,  and  segments 
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FIG.  5  IMPEDANCE  MATCHING  DIAGRAM 
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of  similar  curves,  can  be  used  to  approximate  the  relation  between  P 
and  u  in  a  simple  relief  wave.  These  curves  are  generally  referred 
to  as  'cross  curves'  and  a  curve  can  be  drawn  by  specifying  either  P 
or  u  on  the  Hugoniot,  c.^ve  1. 

In  a  previous  section,  pairs  of  values  of  x  and  Ug  were ^obtained 
from  the  experimental  records.  Because  the  LSGT  calibration^  gives 
u  as  a  function  of  x,  measuring  the  length  of  the  PMMA  attentuator 
yields  the  particle  velocity,  u^,  at  the  PMMA/TNT  interface  just 
before  the  shock  enters  the  TNT.  The  LSGT  calibration  is,  for  a  50/50 
pentolite  donor  (1.56  g/cc) , 

=  1.7735  exp(-0. 01841  x)  t  0.8765  exp(-0.3495  x) 

X  <  36  mm 

(17) 

=  0.0905  +  4.0877  exp(-0. 04451  x) 

X  >  36  mm 

where  ui  is  in  rruti/ysec.  The  pressure  associated  with  u^  is  obtained 
from  equation  (16)  and  is  designated  P]^.  The  particle  velocity  in  the 
TNT,  designated  u^,  is  obtained  from  the  experimental  records  as 
described  in  the  section  "Eleccromagnetic  Velocity  Gage  Techniques". 
There  results  a  set  of  values  of  U]^,  P]^,  and  u^  for  each  experiment. 
It  remains  to  determine  the  value  of  Pg#  the  pressure  actually 
transmitted  into  the  TNT. 

For  each  set  of  data,  the  PMMA  cross  curve  (such  as  curve  2  in 
Figure  5)  is  draw.'  so  that  it  crosses  curve  1  at  u  =  u^.  Then  the 
line  u  =  Ug  is  drawn  and  the  value  of  the  pressure  where  this  line 
intersects  the  cross  curve  is  read.  This  pressure  is  the  value  of  the 
transmitted  pressure,  P^.  This  process  is  repeated  for  each  set  of 
values  of  u^,  P]^,  and  u^. 

Viscoelastic  Effects  in  PMMA 


In  the  previous  section,  the  PMMA  was  assumed  to  behave  in  a 
simple  way  when  it  expanded  from  a  shocked  state.  Rate  effects, 
elastoplastic  and  viscoelastic  effects,  and  the  effects  due  to  the 
increase  of  entropy  behind  the  shock  wave  were  neglected.  Data  on 
the  way  PMMA  behaves  as  it  is  relieved  from  the  shocked  state  are 
scarce.  From  the  work  of  Schuler  and  Nunziato^^  and  Barker  and 
Hollenbach^ ® ,  all  of  Sandia  Laboratories,  it  is  inferred  that  PMMA  is 
a  nonlinear,  viscoelastic  material.  The  data  leading  to  this  infer- 
rence  was  obtained  from  observing  shocked  samples  of  PMMA  with  laser 
interferometry.  These  data  show  that  the  reflected  Hugoniot  as 
discussed  in  the  previous  section  is  only  an  approximation  to  the 
observed  release  path.  In  References  (17)  and  (18) ,  the  observed  data 
were  converted  to  stress  and  strain.  In  the  stress  vs  strain  diagram, 
the  unloading  curve  lies  below  the  loading  (shock  Hugoniot)  curves. 

A  knowledge  of  how  these  experimentally  based  data  appear  in  the 
pressure-particle  velocity  plane  is  needed  for  the  present  work. 
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Unfortunately,  References  (17)  and  (18)  do  not  give  values  of  the 
particle  velocity  for  each  value  of  the  strain  (although  particle 
velocity  was  the  variable  measured  in  the  experiments) . 

The  data  of  References  (17)  '.nd  (18),  however,  can  be  used  to 
obtai.-  values  of  the  particle  velocity  by  using  the  relation 


du  = 


(18) 


which  relates  the  change  in  the  particle  velocity,  u,  to  a  change  in 
the  density,  p,  across  an  increment  of  a  simple  wave.  The  variable 
c  is  the  local  sound  speed.  Engineering  strain,  e,  was  used  in 
References  (17)  and  (18) ,  so  that 

Po 

de  =  _  _  .  (19) 

P  P 

Combining  quations  (18)  and  (19)  yields 


u  =  f£P  de 
J  Po 

The  quantity  c  is  obtained  from 

2  ^  dP  ^  dP  de 

^  3^  3? 


(20) 


(21) 


where  de/dp  is  obtained  from  equation  (19)  to  give 


2  ^  dP  Pq 

c  de  p2  • 

Combining  Equation  (2  0)  and  (22)  gives 

Sets  of  data  from  References  (17)  and  (18)  were  fitted  to  a  cubic  by 
the  method  of  least  squares,  giving  P  as  a  function  of  e,  i.e., 

P  =  A  +  Be  +  Ce2  +  De^  .  (24) 


(22) 


(23) 


Diffe  ‘entiating  Equation  (24)  and  substituting  for  dP/de  in  equation  (23) 
gives  an  integrand  which  can  be  handled  analytically.  There  results 


u 


log[6De  +  2c  R/12D^]|1  +  constant 


(25) 
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where  R  =  (B  +  2Ce  +  3De^)^ .  The  constant  is  evaluated  by  noting  that 
the  material  is  relieved  from  a  state  located  on  the  Hugoniot  for 
which  e  =  and  u  =  ui.  Equation  (25)  is  for  a  forward  facing  wave, 
so  the  cross  curve  is  obtained  from 

u^;,  =  2u;j^  -  u.  (26) 

We  now  have  relations  between  P  and  e  from  the  cubics  fitted  to  the 
data  of  references  (17)  and  (18),  and  bet:»^3en  u^  and  e,  by  equations 
(25)  and  (26) .  Several  of  the  cross  curves  are  plotted  in  Figure  6 
along  with  the  mirror  image  of  the  Hugoniots.  The  numoers  at  the 

upper  ends  of  each  curve  give  the  pressure  in  kilobars  where  the 

curve  intersects  the  Hugoniot  curve  for  forward  facing  shocks.  That 
is,  these  are  the  shock  pressures  from  which  the  PMMA  expands.  For 
8.36,  11.99,  and  18.34  kbar,  there  is  no  significant  difference  between 
the  curves  as  calculated  and  the  reflected  Hugoniots.  For  29.7  kbar, 
the  calculated  curve  lies  below  the  reflected  Hugoniot  by  about  one 
or  two  kilobars  over  the  range  of  pressures  of  8  to  22  kbar.  Hence, 
the  values  of  P^  obtained  by  impedance  matching  should  be  reduced  by 
about  2  kilobars.  Even  though  Schuler  and  Wunziato'^  have  investi¬ 
gated  PMMA  up  to  a  pressure  of  60  kbar,  no  data  above  29.7  kbar  have 

been  found  by  the  authors  with  which  further  calculations  could  be 
made.  Schuler  and  Nunziato^^  do  state  that  viscoelastic  effects 
become  of  secondary  importance  above  "*40  kbar  and  below  ‘^1  kbar. 

In  this  study,  we  have  used  the  calculated  curves  in  Figure  6  to 
estimated  the  correction  due  to  viscoelastic  effects.  Between  30  and 
60  kbars,  the  correction  was  estimated  by  assuming  that  the  visco¬ 
elastic  effect  is  symmetric  about  30  kbar.  Thus  the  correction  becomes 
progressively  smaller  from  30  to  60  kbar. 

Note  that  two  of  the  calculated  curves  in  Figure  6  are  terminated 
before  P  =  0.  This  was  done  because  the  data  of  references  (17)  and 
(18)  also  terminated  at  these  pressures.  Extending  these  curves  would 
require  extrapolating  the  cubic  fit  of  P  vs  e  to  P  =  0. 

Data 


Most  of  the  raw  data  obtained  in  this  study  are  listed  in  Table  1. 
Three  shots  in  the  pressure  range  18-44  kbar  showed  some  evidence  of 
reaction  of  the  explosive  and  are  not  includpd.  The  Pg  data  of  Table  1 
were  obtained  using  the  reflected  PMMA  Hugoniot.  As  is  evident  from 
Table  1,  there  is  some  scatter  in  density  of  the  individual  shots. 

Using  the  computer  code  developed  by  Hermann^  (see  Appendix  A)  to 
predict  the  Hugoniots  of  porous  materials,  the  data  of  Table  1  were 
adjusted  to  the  single  density  of  0.98  g/cc.  This  was  accomplished 
as  follows.  Calculated  Hugoniots  for  TNT  at  0.95,  0.96,  0.97,  0.98, 
and  1.01  g/cc  were  generated  and  plotted.  Then  for  each  experiment, 
the  distance  between  the  calculated  TNT  Hugoniot  at  the  density  of  the 
experiment  and  the  calculated  TNT  Hugoniot  at  0.98  g/cc  along  the 
reflected  PMMA  Hugoniot  was  determined.  The  experimental  datum  point 
was  then  moved  this  distance  along  the  reflected  PMMA  Hugoniot.  This 
correction  for  density  results  in  new  values  for  both  u  and  P  for  the 
TNT  data.  These  results  are  listed  in  Table  2  and  plotted  as  open 
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TABLE  1 

HUGONIOT  DATA  FOR  POROUS  TNT 
(Mirror  Image  Approximation) 


'^e 

mm/]jsec 

Kbar 

Po 

g/cc 

Length 
of  PMMA 
mm 

“i 

mm/ysec 

Pi 

Kbar 

Shot  No. 

0.350® 

3.6 

0.96 

76.2 

0.228 

8.7 

254 

0.450 

8.3 

0.96 

63.5 

0.333 

12.9 

257 

0.440 

8.6 

0.97 

63.5 

0.333 

12.9 

252 

0.700 

12.5 

0.96 

50.8 

0.517 

20.7 

285A 

o.eso*^ 

13.0 

1.01 

50.8 

C.517 

20.7 

251 

0.730 

11.2 

0.97 

50.8 

0.517 

20.7 

253 

0.772"' 

14.2 

0.98 

46.8 

0.560 

22.9 

336 

0.793" 

13.7 

0.98 

47.7 

0.578 

23.8 

333 

0.794" 

13.7 

0.98 

47.7 

0.578 

23.8 

332 

0.970 

17.4 

0.95 

44.5 

0.654 

27.9 

293 

1.080" 

24.5 

0.98 

38.2 

0.837 

38.6 

335 

1.030 

27.8 

0.98 

38.  T 

0.840 

38.8 

290 

1.300 

43.0 

1.01 

25.4 

1.111 

56.9 

316 

1.  310 

42.3 

1.01 

25.4 

1.111 

56.9 

291 

1.340" 

41.0 

0.98 

25.4 

1.111 

56.9 

327 

1.720™ 

55.0 

0.98 

12.7 

1.413 

80.4 

330 

1.720" 

55.0 

0.90 

12.7 

1.413 

80.4 

331 

a)  Peak  of  u,t  curve. 

b)  5  mil  gage,  all  others  1  mil. 
m)  mica  covered  gage  leads. 
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TABLE  2 

HUGONIOT  DATA  FOR  POROUS  TNT  ADJUSTED  TO  0.98  G/CC 
(Mirror  Image  Approximation) 


0.98 

g/cc 

Raw  Data, 

Table  1 

rom/psec 

Pe 

Kbar 

Shot  No. 

Ue 

mm/psec 

Pe 

Kbar 

0.350® 

3.6 

254 

0.350 

3.6 

0.435 

8.4 

257 

0.450 

8.3 

0.440 

8.6 

252 

0.440 

8.6 

0.680 

13.0 

285A 

0.700 

12.5 

0.700^ 

12.5 

251 

0.680 

13.0 

0.725 

11.4 

253 

0.730 

11.2 

0.772™ 

14.2 

336 

0.772 

14.2 

0.793™ 

13.7 

333 

0.793 

13.7 

0.794™ 

13.7 

332 

0.794 

13.7 

0.960 

18.0 

293 

0.970 

17.4 

1.080™ 

24.5 

335 

1.080 

24.5 

1.030 

27.8 

290 

1.030 

27.8 

1.330 

41.5 

316 

1.300 

43.0 

1.340 

41.0 

291 

1.310 

42.3 

1.340™ 

41.0 

327 

1.340 

41.0 

1.720™ 

55.0 

330 

1.720 

55.0 

1.720™ 

55.0 

331 

1.720 

55.0 

a)  Peak  of  u,t  curve, 

b)  5  mil  gage,  all  others  1  mil. 
m)  Mica  covered  gage  leads. 
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circles  in  Figure  7.  The  changes  are  small  except  for  the  points  at 
-  40  kbar  where  the  correction  *.s  the  largest. 

The  data  contain  five  sets  of  replicate  shots  (i.e.,  each  set  has 
the  same  attenuator  length).  For  the  first  set,  shots  252  and  257,  u 
has  the  values  of  0.440  and  0.435  mm/ysec  (see  Table  2).  The  resultant 
values  of  Pg,  8.6  and  8.4  kbar,  are  within  ±  1.2%  of  the  average  for 
this  attenuator  length.  For  the  second  set,  shots  285A,  251,  and  253, 
Ug  =  0.680,  0.700,  and  0.725  mm/ysec,  respectively,  with  resulting 
P^*s  of  13.0,  12.5,  and  11.4  kbar.  These  pressures  are  within  ±  7.3% 
of  the  average  at  this  attenuator  length.  For  the  third  set,  shots 
332  and  333,  Ug  =  0.794  and  0.793  ram/ysec,  respectively,  with  the  same 
value  of  Pg,  13.7  kbar.  For  the  fourth  set,  shots  316,  291,  and  327, 
Ug  =  1.330,  1.340,  and  1.340  mm/ysec,  respectively,  with  resulting 
Pg*s  of  41.5,  41.0,  and  41.0  kbar.  The  pressures  are  within  ±  0.8% 
of  the  average.  The  final  set  of  replicate  shots,  330  and  331,  have 
the  same  value  of  Ug,  1.720  mm/ysec,  and  Pg,  55.0  kbar. 

Because  the  calculated  release  paths  in  Figure  6  do  not  coincide 
with  the  reflected  Hugoniots,  we  have  to  repeat  the  impedance  matching 
solutions  which  were  described  previously.  This  time  we  use  infor¬ 
mation  such  as  that  shown  in  Figure  6  to  draw  viscoelastic  release  paths 
through  each  predetermined  point,  (ui,Pj^)  ,  i.e.,  points  such  as  that 
labeled  A  in  Figure  6.  Each  of  these  points  lie,  of  course,  on  the 
Hugoniot  for  forward  facing  shocks  in  the.  PMMA  and  are  determined 
solely  by  the  length  of  the  attenuator,  see  Equations  16  and  17.  The 
amount  of  the  offset  of  the  new  curves  from  the  reflected  Hugoniots 
is  determined  by  interpolating  in  the  data  of  Figure  6  and  assuming 
that  the  correction  is  symiTietrical  about  the  release  path  for  a  30  kbar 
shock.  This  process  gives  corrections  which  are  significant  for 
values  of  P^  which  lie  between  10  and  40  kbars,  see  Table  3.  The 
maximum  correction  in  Pg,  1.5  kbars,  occurs  at  a  value  of  Ug  of  about 
1.05  mm/ysec  and  Pg  of  about  25  kbars.  On  a  percentage  basis,  the 
maximum  correction,  about  9%,  occurs  at  0.725  mm/ysec  and  11.4  kbars, 
see  Table  3. 

Data  on  the  dynamic  compressibility  of  1.0  g/cc  TNT  has  been 
published  by  two  groups  of  Russian  investigators.  Vasil 'ev, 
Bolkhovitinov,  and  Khrisotoforov^® employed  the  EMV  gage  in  1.0  g/cc 
TNT.  The  P,u  data  obtained  from  their  results  are  plotted  in  Figure  8 
as  open  triangles.  The  results  of  Table  3  are  plotted  in  Figure  8  as 
solid  circles  for  comparison.  The  errors  in  pressure  of  our  data 
were  computed  using  the  method  described  in  Appendix  B.  The  error  in 
particle  velocity  lies  within  the  size  of  the  circles.  As  is  evident 
from  the  figure,  Vasil *ev*s  results  diverge  from  ours  as  the  pressure 
increases.  This  divergence  is  due  to  the  fact  that  Vasil’ev's  EMV 
gage  is  located  inside  the  TNT  sample  and  is  thus  affected  by  any 
reaction  which  has  had  time  to  build  up  in  the  TNT  as  the  shock  wave 
travels  from  the  interface  to  the  gage  location.  These  authors  are 
aware  of  this  because  they  report  that  the  shock  velocity  builds  up 
with  increasing  distance  in  the  sample.  Thus,  except  for  the  very  low 
pressure  points,  their  data  represents  points  on  various  reacted 
Hugoniots.  Dremin,  Koldunov,  and  Shvedov^®  have  investigated  the 
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TABLE  3 

HUGONIOT  DATA  FOR  POROUS  TNT,  Pq  =  0.98  G/CC, 
ADJUSTED  FOR  ESTIMATED  VISCOELASTIC  EFFECTS  IN  PMMA 


«e 

Pe 

U 

V 

mm/ysec 

Kbar 

mm/ysec 

cc/g 

Shot  No 

0.350® 

3.6 

1.05 

0.680 

254 

0.435 

8.4 

1.97 

0.795 

257 

0.440 

8.6 

1.99 

0.795 

252 

0.680 

12.0 

1.80 

0.635 

z  b 

0.700^ 

11.5 

1.68 

0.594 

251 

0.725 

10.4 

1.46 

0.515 

253 

0.772"* 

13.2 

1.74 

0.551 

336 

0.793™ 

12.7 

1.63 

0.525 

333 

0.794™ 

12.7 

1.63 

0.525 

332 

0.960 

17.0 

1.80 

0.478 

293 

1.080™ 

23.0 

2.17 

0.513 

335 

1.030 

26.3 

2.60 

0.617 

290 

1.330 

40.5 

3.11 

0.583 

316 

1.340 

40.0 

3.06 

0.572 

291 

1.340™ 

40.0 

3.05 

0.572 

327 

1.720" 

55.0 

3.26 

0.428 

330 

1.720"* 

55.0 

3.26 

0.428 

331 

a)  Peak  of  u,t  curve. 

b)  5  mil  gage,  all  others  1  mil. 
ra)  Mica  covered  gage  leads . 
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—  EQ.  (9) 


PARTICLE  VELOCITY,u,  mm/jjsec 


FIG.  8  COMPARISON  OF  DATA  FOR  POROUS  TNT 
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buildup  of  particle  velocity  and  shock  velocity  in  1.0  g/cc  TNT  using 
the  EMV  gage.  Although  most  of  their  data  were  obtained  from  gages 
located  in  the  TNT  charges,  Dremin  et  al  report  some  results  where  the 
gage  was  located  at  the  interface.  Such  data  will  lie  on  the  unreacted 
Hugoniot  as  was  discussed  in  the  earlier  section  entitled  "Effects  of 
Reaction  on  Gage  Output”.  These  interface  data,  which  were  obtained 
from  Figure  2  of  Dremin *s  paper,  are  plotted  in  Figure  8  as  open 
squares.  As  is  evident  from  the  figure,  Dremin *s  data  agree  very  well 
with  ours. 

A  Hugoniot  is  obtained  experimentally  by  measuring  two  of  the 
following  quantities:  pressure,  particle  velocity,  shock  velocity, 
and  volume.  Once  two  of  these  quantities  have  been  related  to  each 
other,  the  other  two  can  be  calculated  using  hydrodynamic  relations. 

In  this  study,  pressure  and  particle  velocity  were  obtained  experi¬ 
mentally.  Besides  the  pressure  -  particle  velocity  plane,  Hugoniots 
are  frequently  presented  in  the  pressure  -  volume  and  shock  velocity  - 
particle  velocity  planes. 

T'.ie  shock  velocity  (U)  -  particle  velocity  (u)  representation  of 
the  Hugoniot  for  0.98  g/cc  TNT  is  obtained  by  using  u«  and  from 
Table  3  and  Equation  (6)  with  =  1/Po  ^o  assumed  to  be® zero. 

The  resulting  shock  velocity  values  are  presented  in  column  3  of 
Table  3  and  plotted  as  a  function  of  u  in  Figure  9. 

The  pressure  -  volume  (P-v)  representation  of  the  Hugoniot  for 
0.98  g/cc  TNT  is  obtained  by  using  the  Pg  and  u^  data  of  Table  3  and 
the  following  equation: 


V 


^of  - 


(27) 


where  v^f  is  the  initial  volume  of  the  TNT.  The  calculated  values  of 
the  volume  are  presented  in  column  4  of  Table  3.  Figure  10  is  the 
resulting  P,v  plot.  The  errors  in  P  and  v  were  obtained  using  Equa¬ 
tions  (B7)  and  (B8)  of  Appendix  B.  The  large  errors  in  volume  is  due 
to  the  fact  that  even  small  errors  in  P  and  u  can  cause  a  relatively 
large  error  in  v. 

COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  HUGONIOTS  FOR  POROUS 
EXPLOSIVES,  DATA  FROM  NAVAL  SURFACE  WEAPONS  CENTER  AND  OTHER  SOURCES 

The  experimentally  determined  Hugoniot  for  TNT  at  0.98  g/cc  can  now 
be  compared  with  the  results  from  computational  models.  The  required 
data  for  Herrmann’s  code  are:  the  Hugoniot  of  the  solid  material,  a 
value  of  r (v) (or  Tq  if  f (v)  is  not  available),  and  the  specific  volume 
of  the  solid.  The  crushup  part  of  the  code  cannot  be  used  because 
values  of  important  parameters  are  not  known,  i.e.,  the  elastic  limit, 
the  crushup  pressure,  and  the  sound  speed.  The  Hugoniot  and  Vq  for 
solid  TNT  at  a  density  of  1.614  g/cc  are^^ 

U  =  2.39  +  2.05U  (28) 

and 

To  =  0.737  (29) 

U  and  u  are  in  mm/ysec.  The  resulting  Hugoniot  for  unreacted 


where 
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FIG.  9  U,u  PLOT  Or  0.98  g/cc  TNT  HUGONIOT 


25 


UNCLASSIFIED 


•  THIS  WORK 


UNCLASSIFIED 

NOLTR  74-213 

TNT  at  a  density  of  0.98  g/cc  is  shown  as  the  solid  curve  in  Figure  8 
(P,u  plane).  Except  for  the  two  pressure  points  at  -  8.5  kbar,  the 
calculated  Hugoniot  agrees  with  the  experimental  data  to  within  the 
calculated  error  in  pressure  and  particle  velocity.  The  data  at 
-8,5  kbar  are  better  explained  in  the  P-v  plane  and  will  be  discussed 
below.  Another  computed  Hugoniot  is  shown  in  Figure  8  as  a  dashed 
line.  This  is  based  on  the  theory  given  in  Reference  (9)  and  is 
expressed  as  Equation  (9)  above.  It  requires  the  use  of  Equations  (6) 
and  (28).  The  results  of  this  model  do  not  represent  the  experimental  data 
at  all  well.  We  have  not  tried  to  readjust  the  values  of  any  of  the 
parameters  in  Equation  (9)  to  get  a  better  fit  to  the  experimental 
data.  The  Herrmann  model  gives  such  good  results  with  no  adjustment 
of  parameters  that  only  a  comparison  of  this  model  and  the  experimental 
results  in  the  U,u  and  P,v  planes  will  be  presented  in  the  following 
discussion. 

The  calculated  Hugoniot  for  0.98  g/cc  TNT  in  the  U,u  plane  is 
plotted  in  Figure  9  as  a  solid  curve.  Except  for  the  two  data  points 
at  u  -  0.44  mm/psec  and  U  -  1,98  mm/psec,  the  calculated  Hugoniot 
agrees  with  the  experimental  data.  As  mentioned  above,  these  two 
data  points  are  more  profitably  discussed  in  the  P,v  plane.  The 
experimental  data  (omitting  the  above  mentioned  two  data  points)  were 
fitted  to  a  straight  line  by  the  method  of  least  squares  and  the 
resulting  equation  is 

U  =  0.366  +  1.813  u  (30) 

where  U  and  u  are  in  mm/psec.  The  calculated  Hugoniot  for  this  density 
is  nearly  a  straight  line.  The  constants  for  this  straight  line  were 
determined  by  fitting  the  computed  results  to  a  straight  line  by  the 
method  of  least  squares.  (The  computed  results  are  individual  points. 
The  number  of  points  and  the  separation  of  the  points  is  determined  by 
the  step  size  used  in  the  computer  code,  see  Appendix  A.)  The  resulting 
equation  is 

U  =  0.307  +  1.848  u.  (31) 

The  experimental  result  is  essentially  parallel  to  but  lies  0.059 
mm/psec  above  the  calculated  Hugoniot. 

In  the  P,v  plane  (Figure  10) ,  the  calculated  Hugoniot  is  again 
plotted  as  a  solid  curve.  The  agreement  between  the  experimental  data 
and  the  calculated  Hugoniot  does  not  seem  to  be  as  good  as  in  the  P,u 
and  U,u  planes.  But  as  was  mentioned  in  the  preceeding  section,  a 
large  error  in  the  value  of  v  can  result  from  small  errors  in  P  and  u. 
Taking  this  fact  into  account,  the  agreement  is  reasonable. 

The  two  data  points  at  P  ~  8.5  kbar  and  v  =  0.795  cc/g  of  Figure  10 
correspond  to  the  two  data  points  in  the  P,u  and  U,u  planes  which  did 
not  agree  with  the  calculated  Hugoniot  (Figures  8  and  9) .  The  calcu¬ 
lated  Hugoniot,  however,  shown  in  Figures  8,  9,  and  10  was  obtained 
by  neglecting,  among  other  phenomena,  the  crushup  process.  These  two 
data  points  probably  lie  on  the  crushup  portion  of  the  Hugoniot.  The 
dashed  line  in  Figure  10  represents  a  possible  crushup  path  of  the 
Hugoniot.  It  is  an  arbitrary  straight  line  originating  at  Vof  and 
ending  at  the  calculated  Hugoniot.  However,  there  are  too  few  data  • 
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available  to  properly  define  this  portion  of  the  Hugoniot.  Another 
possibility  which  could  explain  the  position  of  these  two  data  points 
is  that  melting  of  the  TNT  may  be  occurring.  This  possibility  has  Leen 
investigated  but  no  conclusion  about  melting  could  be  drawn  because 
buildup  to  detonation  occurred  too  rapidly  and  erased  any  hint  of  a 
phase  change.  The  subject  of  melting  of  the  TNT  is  more  full\  discussed 
in  Appendix  C.  The  investigation  of  melting  involved  experimental  work 
with  1.3  g/cc  TNT  which  yielded  some  data  for  the  unreacted  state. 

These  results  are  presented  in  Appendix  D.  A  third  possibility  for  the 
position  of  these  two  data  points  is  that  crush-up  and  melting  could 
be  occurring  simultaneously.* 

Data  on  five  other  porous  explosives  are  available  in  the  Russian 
literature  from  which  the  unreacted  Hugoniots  can  be  calculated  and 
compared  with  experimental  results.  Densities  of  the  porous  explosives, 
the  coefficients  for  the  relation  U  =  Cq  +  Su  for  the  solid  materials 
and  values  of  Vq  for  each  explosive  are  given  in  Table  4.  Tabulated 
results  of  the  computations  using  the  Herrmann  model  are  listed  in 
Table  A4  of  Appendix  A.  Some  of  these  results  are  discussed  below. 

Data  for  NB-40  (60%  pyroxyline,  40%  nitroglycerine)  are  from 
Reference  22;  the  remainder  are  from  Reference  10.  Computed  and  exper¬ 
imental  results  are  shown  for  NB-40  in  Figure  11  as  a  plot  of  shock 
velocity  as  a  function  of  particle  velocity.  The  experimental  results 
are  more  closely  approximated  by  the  results  from  Herrmann's  code 
(solid  curve)  than  by  those  from  Equation  (9)  (the  Russian  approximation, 
dashed  curve) .  This  makes  two  cases  in  which  this  Russian  approxima¬ 
tion  gives  the  worse  set  of  results.  In  the  following,  only  the 
Herrmann  model  will  be  used  because  of  its  demonstrated  superiority. 

For  the  other  four  materials  listed  in  Table  4,  results  are  shown 
in  the  U,u  plots  of  Figures  12  through  15.  Results  shown  are  (a)  the 
experimental  results  as  reported  in  Reference  10,  (b)  results  of 
computations  as  given  by  the  same  source  using  the  values  of  Table  4, 
and  (c)  results  computed  with  the  Herrmann  model  also  using  the  values 
of  the  parameters  given  in  Table  4.  In  their  computations,  the  Russian^ 
used  the  method  based  on  the  0®K  isotherm  as  briefly  described  in  the 
section,  "Computational  Models".  The  Russians  reported  computed 
results  at  25,  50,  75,  100,  and  150  kbar.  These  data  are  recorded  in 
Figures  12  through  15  ar  points  rather  than  as  curves. 

The  results  for  RDX  are  shown  in  Figure  12  where  it  is  seen  that 
the  Russian  computations  are  in  no  better  agreement  with  their  experi¬ 
mental  data  than  are  the  results  of  Herrmann’s  model.  Thus,  it  appears 
that  the  much  more  involved  method  used  by  the  Russians  is  not  neces¬ 
sary  because  the  simpler  model  gives  equivalent  results. 

For  tetryl,  PETN,  and  AN,  the  Russians  list  two  values  of  the  shock 
velocity  for  each  listed  particle  velocity.  They  did  this  because  shock 
velocity  was  determined  by  two  different  methods.  Their  directly 
measured  values  of  shock  velocity  for  these  three  materials  have  been 
plotted.  Where  possible,  the  value  which  they  obtained  by 

♦Data  on  the  buiL^up  oT  pirtlcXe  velocity  in  0.98  and  1.30  g/cc  TNT 
are  presented  in  Appendix  E. 
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TABLE  4 

PARAMETERS  FOR  FIVE  EXPLOSIVES  FOR  THE  HERRMANN  MODEL 


Density  Density  of  Solid 


Explosive 

q/cc 

mm/ysec 

S 

g/cc 

_  Re^ 

NB-40 

1.00 

1.70 

1.85 

1.50 

1.63 

22 

RDX 

1.00 

2.87 

1.61 

2.60 

1.80 

10 

Tetryl 

0.86 

2.17 

1.91 

1.65 

1.73 

10 

PETN 

0.82 

2.42 

1.91 

1.70 

1.77 

10 

AN 

0.86 

2.20 

1.96 

0.90 

1.73 

10 
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what  they  c^ll  the  "reflection  technique"  (i.e.,  impedance  matching) 
is  also  plotted.  They  explain  how  the  coincidence  of  the  results  of 
these  two  methods  assures  them  that  reaction  has  not  affected  their 
results.  For  tetryl,  Figure  13,  both  computational  models  fit  the 
experimental  data  very  well.  For  PETN,  Figure  14,  our  computations 
give  values  of  U  which  are  too  small  for  given  values  of  u,  while  the 
Russian  computations  follow  their  data.  Again,  we  have  not  tried 
adjusting  parameters  to  make  the  Herrmann  model  fit  the  da. a  more 
closely.  For  AN,  Figure  15,  both  computations  give  values  of  U  which 
are  too  small. 

The  above  comparisons  are  summarized  in  Table  5.  The  materials 
are  listed  in  the  left  hand  column  of  the  table.  The  next  column 
gives  the  density  of  the  porous  material  used  in  the  experimental  work. 
The  next  three  columns  indicate  (1)  how  our  computations  agree  with  the 
experimental  data,  (2)  how  the  Russian  calculations^® compare  with  their 
own  experimental  data,  and  (3)  how  our  computations  compare  with  the 
Russian  computations. 


TABLE  5 

COMPARISON  OF  COMPUTED  AND  EXPERIMENTAL  RESULTS  FOR  SIX  EXPLOSIVES 

Density  Our  Calc.  vs.  Russian  Calc.  vs.  Our  Calc.  vs. 
Material  g/cc  Experimental  Data  Experimental  Data  Russian  Calc. 

TNT  0.98  good  +  + 

NB*  1.00  good  +  + 

RDX*  1.00  good  good  good 

tetryl*  0.86  good  good  good 

PETN*  0.82  poor  good  poor 

AN*  0.86  poor  poor  good 

*  Experimental  work  from  Reference  10. 

+  No  results  reported  in  Russian  literature. 


We  conclude  that  our  calculations  fit  the  data  at  least  as  well 
as  do  the  Russian  calculations,  which  are  much  more  involved.  Hence, 
we  see  no  need  to  use  the  more  complicated  computational  method  used 
by  the  Russians.  For  AN,  the  computations,  in  both  cases,  may  be  in 
error  because  of  errors  in  the  Hugoniot  of  the  solid  AN.  This  material 
changes  phases  under  pressure  so  that  the  Hugoniot  cannot  be  represented 
by  a  single  linear  relation  between  shock  and  particle  velocities.  Of 
course,  the  experimental  data  for  AN,  and  for  PETN,  may  also  be  in 
error  even  though  the  Russians  tried  to  avoid  letting  reaction  affect 
their  results. 
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Up  to  this  point  in  this  report,  we  have  used  the  Herrmann  code 
to  calculate  the  0.98  g/cc  Hugoniot  for  TNT  only.  Figures  16  and  17 
give  computed  Hugoniots  for  TNT  for  densities  from  1.0  to  1.60  g/cc.* 

In  Appendix  D  we  compare  experimental  results  with  the  computed 
Hugoniot  for  TNT  at  1.3  g/cc.  Here  we  note  that  in  Figure  16  the 
relation  for  1.0  g/cc  TNT  is  very  nearly  linear  —  this  is  why 
equation  (31)  can  be  used.  At  greater  initial  densities  (smaller 
values  of  a©),  the  relations  tend  to  become  curved  for  small  values  of 
the  particle  velocity,  that  is,  at  low  pressures.  This  is  noticeable 
even  for  1.60  g/cc  TNT.  Because  we  have  not  been  able  to  compute  the 
shock  variables  in  the  crushup  regime,  we  cannot  attach  a  great  deal 
of  importance  to  the  curvature  mentioned  above.  That  is,  if  we  had 
been  able  to  compute  the  values  of  these  variables  in  the  crushup 
regime,  the  curvature  may  have  been  more  or  less  than  we  show  in 
Figures  16  and  17.  Actually,  our  U,u  curves  resemble  those  published 
by  Herrmann^  in  his  Figure  2  which  is  for  aluminum.  In  his  calcula¬ 
tions,  he  neglected  the  elastic  response  of  the  porous  aluminum.  For 
iron,  the  curves  become  more  complex,  see  Figure  6  in  Reference  (7) . 

In  that  work,  Herrmann  tried  to  account  for  the  elastic  response  of 
the  distended  iron,  and  the  crushup  process.  Hence,  we  should  expect 
to  obtain  curves  that  are  more  complicated  than  those  shown  in  our 
Figure  16  when  we  are  able  to  account  for  the  entire  compaction  process. 


CONCLUSIONS  AND  FUTURE  WORK 

Our  results  show  that  Hugoniots  for  porous  TNT  can  be  calculated 
to  a  good  degree  of  accuracy  by  using  the  mathematical  model  described 
by  Herrmann’.  Most  of  our  work  was  with  TNT  at  0.98  g/cc;  a  small 
amount  with  TNT  at  1.3  g/cc.  Furthermore,  we  reproduce  experimental 
Hugoniots  given  in  the  Russian  literature^®'’’  for  NB.  RDX,  and  tetryl 
to  a  good  degree  of  accuracy.  For  PETN  and  for  AN,  our  computations 
agree  less  well  with  the  Russian  data^®'”.  This  failure  may  be  due 
to  inaccuracies  in  the  Russian  Hugoniot  data  on  the  low  density  forms 
of  each  material,  or  on  the  values  of  the  parameters  which  the  Russians 
use  to  describe  the  Hugoniots  of  the  solids  or  a  combination  of  these. 
The  Russians  themselves  do  not  obtain  close  agreement  between  calculated 
and  experimental  Hugoniots  for  AN.  This  material  is  hygroscopic  which 
makes  it  difficult  to  keep  pure  in  experimental  work.  Also,  it  under¬ 
goes  phase  changes  under  pressure.  For  these  reasons,  the  Hugoniots 
of  porous  AN  may  be  difficult  to  model.  In  the  case  of  PETN,  a  fairly 
sensitive  explosive,  reaction  may  have  caused  errors  in  the  Russian 
experiments,  especially  with  the  porous  forms  of  the  explosive.  These 
errors  may  exist  in  spite  of  the  prec;autions  employed  by  the  Russians 
(see  the  previous  section  of  this  report) . 


*Tables  A2  and  A3  in  Appendix  A  give  values  of  some  of  the  flow 
variables  as  functions  of  the  specific  volume,  v,  for  TNT  at  densities 
ranging  from  1.0  g/cc  to  1.614  g/cc. 
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FIG.  14  COMPARISON  OF  EXPERIMENTAL  AND  COMPUTED  HUGONIOTS  FOR  PETN  AT  0.82  g/cc 
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V  SHOCK  VELOCITY  MEASURED  DIRECTLY,  REF.  10 
A  SHOCK  VELOCITY  OBTAINED  GRAPHICALLY,  REF.  10 
O  COMPUTED  HUGONIOT,  REF.  10 
- HERRMANN'S  CODE,  REF.  7 


PARTICLE  VELOCITY,u,  mrr/Msec 


FIG.  15  COMPARISON  OF  EXPERIMENTAL  AND  COMPUTED  HUGONIOTS  FOR  AN  AT  0.86  g/ 


36 

UNCLASSIFIED 


ec 


UNCLASSIFIED 
NOLTR  74-213 


FIG.  17  PRESSURE  VS  PARTICLE  VELCXITY  FOR  TNT  AT  DIFFERENT  DENSITIES 
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The  successes  which  we  have  reported  above  lead  us  to  conclude 
that  the  mathematical  model  which  we  are  using  is  adequate  for 
predicting  good  approximations  for  the  Hugoniots  of  porous  explosives. 

It  would  be  most  surprising  if  the  outcome  had  been  negative.  The 
model,  and  variations  of  the  model,  have  been  used  successfully  for 
years  to  describe  the  compaction  of  porous,  unreacting  solids  by  shock 
waves.  The  difficulty  in  working  with  explosives  is,  of  course,  their 
tendency  to  react  when  shocked.  This  challenges  the  experimentalist 
when  he  tries  to  acquire  Hugoniot  data  for  explosives  in  either  their 
solid  or  porous  form. 

Otherwise,  the  degree  of  success  in  computing  the  Hugoniot  for  a 
porous  explosive  depends  on  the  same  factors  as  for  an  inert  material. 
For  example,  if  the  explosive  undergoes  a  phase  change  under  pressure 
(hydrostatic  or  by  shock) ,  the  model  cannot  be  applied  in  the  form  now 
available  to  us.  For  example,  the  low  melting  point  of  TNT  caused 
some  ambiguity  in  the  present  study.  If  the  solid  material  does  not 
undergo  a  phase  change  under  pressure  (and  the  measurements  are  made 
before  reaction  can  proceed  to  any  significant  degree  in  the  experiment) 
the  relation  between  the  shock  and  particle  velocities  is  apt  to  be 
linear  as  required  by  the  simplified  computer  code  reported  in 
Appendix  A.  Curvature  in  the  U,u  relation  for  the  solid  can  be  handled 
by  using  a  power  series  which  relates  the  pressure  to  the  volumetric 
strain,  n  =  (vq  -  v)/Vq,  see  reference  (7). 

The  Gruneisen  parameter  of  the  solid  material  must  be  known. 

Studies  of  the  same  material  in  both  solid  and  porous  forms  gives,  in 
principle,  values  of  the  Gruneisen  parameter,  F.  Our  TNT  data  are  too 
imprecise  and  too  restricted  in  range  for  us  to  get  any  improvement  in 
the  value  of  F  for  TNT,  The  Russians,  using  this  method,  have  given 
values  of  F  for  some  explosives  to  pressures  as  great  as  200  kbar. 

This  is  an  extraordinary  feat  experimentally  in  view  of  the  rapidity 
of  the  onset  of  reaction  in  most  of  these  materials  at  these  higher 
pressures.  The  value  of  Fq  can  be  inferred  from: 

To  =  ac^/Cp 

wliare  a  is  the  volumetric  thermal  expansion  coefficient,  c  is  the 
ambient  sound  speed  and  Cp  is  the  specific  heat  at  constant  pressure. 

For  many  solid  explosives,  reliable  values  of  these  parameters  do  not 
exist.  Hence,  determinations  of  Cp,  c^,  and  a  are  necessary  supple¬ 
mentary  work  for  successful  use  of^the  Herrmann  Code. 

A  big  part  of  the  problem  in  the  study  of  shock  compression  of 
porous  media  is  what  takes  place  when  low  amplitude  shocks  propagate 
through  the  material.  We  have  not  been  able  to  deal  with  this  problem 
in  our  TNT  work,  partly  because  of  the  experimental  limitation  which 
prevents  us  from  producing  low  pressure  shocks  accurately  and  reproduc- 
ibly.  The  potential  simultaneous  melting  of  the  TNT  also  contributes 
to  the  problem  of  interpreting  and  using  the  data.  Application  of 
the  mathematical  models  to  the  compaction  region  must  wait  until  data 
are  available  on  the  elastic  behavior  of  the  porous  materials.  (Some 
of  the  relevant  data  can  be  obtained  from  sonic  experiments.)  Another 
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important  parameter  is  the  pressure  at  which  the  voids  are  essentially 
completely  closed.  For  TNT,  this  measurement  may  hcive  been  obscured, 
we  believe,  by  melting  occurring  simultaneously  with  the  completion  of 
crushup.  It  is  most  unfortunate  that  we  do  not  have  the  data  necessary 
to  apply  the  computer  models  to  the  crushup  process  in  explosives 
because  the  crushup  pressure  may  be  about  the  same  as  the  critical 
pressure  of  the  explosive  in  the  large  scale  gap  test  (LSGT)  for  the 
very  porous  charges.  Hence,  we  cannot  convert  our  values  of  Pg  (the 
pressure  in  the  PMMA  attenuator)  to  values  of  Pi,  the  critical  pressure 
in  the  porous  explosive  because  we  could  not  obtain  a  reliable  non¬ 
reactive  explosive  Hugoniot  in  this  low  pressure  region. 

No  future  work  on  the  non-reacted  Hugoniot  of  porous  explosives 
is  planned  for  two  reasons.  The  first  reason  is  that  we  have  esta¬ 
blished  the  applicability  of  Herrmann’s  code  to  within  the  precision 
of  our  own  (and  the  Russian)  measurements  in  the  pressure  region  of 
10  -  55  kbar  in  this  work  and  10  -100  kbar  in  the  Russian  work.  The 
second  reason  is  that  the  present  experimental  procedures  are  inade¬ 
quate  for  studying  the  crushup  region,  i.e.,  the  region  below  10  kbar. 
The  Russian  researchers  may  also  have  problems  in  this  region  because 
they  do  not  publish  any  data  below  10  kbar. 

Future  work  will  involve  the  determination  of  non-reactive 
Hugoniots  of  voidless  or  near  voidless  high  explosives  for  which  no 
such  data  are  now  available.  The  explosives  which  will  be  studied 
should  satisfy  one  or  more  of  the  following  conditions. 

1)  The  explosive  dead-presses.  That  is,  it  will  not  propagate 
a  detonation  in  the  present  experimental  configuration, 

2)  A  value  for  Tq  (or  the  function  T (v) )  is  available  or  can  be 
determined  for  the  explosive. 

3)  The  explosive  is  used  in  the  physical  modeling  of  existing 
explosive  mixtures  or  propellants. 

SUMMARY 

A  theory  of  the  response  of  porous  materials  to  shock^'^'^  which 
was  found  to  apply  to  non-reactive  materials,  has  been  extended  to  a 
potentially  reactive  material,  porous  TNT.  A  computer  code  developed 
by  Herrmann^  and  based  on  this  theory  used  to  Calculate  Hugoniots  for 
unreacted  TNT  for  the  range  of  densities  from  1.60  to  0.98  g/cc.  The 
calculations  were  confirmed  by  obtaining  an  experimental  Hugoniot  for 
porous  TNT  at  a  density  of  0.98  g/cc.  The  calculated  and  experimental 
Hugoniots  agree  to  within  2  kbar  over  the  range  0.35-1.72  mm/ysec  and 
3.6-55.0  kbar. 

Thus,  it  appears  that  Hugoniots  for  unreacted,  porous  explosives 
can  be  obtained  using  a  computational  method  (provided  that  a  good 
Hugoniot  is  available  for  the  near  voidless  unreacted  explosive)  with 
a  high  degree  of  confidence.  This  method  results  in  a  saving  of  time 
and  money  compared  to  the  uso  of  experimental  method  to  obtain  these 
Hugoniots. 
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APPENDIX  A 

COMPUTER  CODE  FOR  CALCULATING  HUGONIOTS  FOR  POROUS  SOLIDS 

The  computer  code  given  by  Herrmann  has  the  option  for  computing 
the  pressure  of  a  porous  material  in  the  crushup  region.  Actually, 
the  major  part  of  that  code  consists  of  the  crushup  option.  All 
that  is  needed  for  the  work  which  we  have  reported  is  a  means  of 
computing  values  from  Equations  4  and  8  of  the  main  text.  A  code  for 
doing  this  is  given  in  Table  Al.  It  is  written  in  CDC  Fortran,  ard 
is  designed  to  run  from  a  teletype  terminal  using  CDC  "INTERCOM" 

System.  Input  to  the  code  is  accomplished  by  using  an  "INPUT"  state¬ 
ment.  The  code  asks  for  input  data  as  shown  on  the  first  line  in 
Table  A2.  These  are  values  for  Cq  and  S,  see  Equation  5,  ,  the 

initial  porosity  parameter,  Uq/  i.e.,  VQf/vQg,  and  finally  Vqq.  In 
the  code,  these  variables  are  named  CO.  S,  GO,  AO,  and  VS,  respectively. 
The  data  are  immediately  printed  out  by  the  computer  for  checking  by 
the  operator.  Three  quantities  are  given  on  the  next  line  in  Table  2, 
The  first  is  the  value  of  VS,  the  specific  volume  of  the  solid.  The 
second  is  that  at  which  Pf  goes  to  infinity,  see  equation  4  in  the 
main  text.  The  difference,  VS  -  V2,  is  divided  by  40  to  give  an 
increment  in  the  specific  volume  for  the  following  computation.  This 
increment,  V3,  is  used  in  a  "DO"  loop  to  increment  the  specific  volume 
20  times,  or  until  PP,  i.e.,  Pf,  exceeds  0.1  megabar  (for  some  runs, 
the  test  value  was  0.175  megabars). 

The  tabular  output  of  the  code  consists  of: 

V,  the  specific  volume  of  the  shocked  material,  i.e.,  v, 

PH,  the  pressure  in  the  solid  at  specific  volume,  i.e.,  Psr 

PP,  the  pressure  in  the  porous  material  at  specific  volume  v, 
i.e.,  Pf, 

U8,  the  shock  velocity,  i.e.,  U, 

U9,  the  particle  velocity,  i.e.,  u, 

El,  the  change  in  the  internal  energy  due  to  the  shock,  given  by 
Pf(Vof  -  V)x  1.195  X  lO^. 

The  change  in  internal  energy  (El  above)  can  be  used  to  obtain  a  first 
approximation  to  the  temperature  jump  across  the  shock  by  dividing  by 
the  specific  heat.  When  a  table  is  completed,  the  computer  asks  for  a 
new  set  of  data. 

Units  used  in  the  code  are  cm/usec  for  velocities  (in  contrast  to 
the  use  of  mm/psec  in  the  main  text) ,  cm^/g  for  volume,  megabars  for 
pressure  and  cal/g  for  energy. 
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Table  A3  gives  the  results  of  the  code  for  TNT  at  densities  of 
1.07,  1.25,  1.42,  1.49,  1.58,  1.60,  and  1.614  g/cc.  The  last  density 
is  that  of  TNT  in  what  we  regard  in  this  report  as  the  solid  form. 
Consequently,  the  values  of  PH  and  PP  (i.e.,  Pg  and  Pf)  are  identical 
for  the  same  value  of  v.  Part  of  the  results  given  in  Table  A3  are 
displayed  in  Figures  16  and  17  in  the  main  text. 

Table  A4  gives  the  results  of  the  computations  on  which  Figures  11 
through  15  are  based;  the  input  data  and  its  sources  are  given  in 
Table  4  of  the  main  text.  An  interesting  feature  of  Table  A4  is  the 
results  for  RDX^at  1.0  g/cc.  The  increment  in  the  volume,  Av,  turned 
out  to  be  negative  -  the  code  extrapolates  the  solid  Hugoniot  into  the 
tension  region,  hence  the  values  of  PH  are  negative.  The  values  of 
PP  (i.e.,  Pf)  are  positive,  so  that  the  v  vs  Pf  curve  has  positive 
slope  rather  than  negative  slope.  This  behavior  is  inherent  in  the 
model  and  has  been  observed  r  perimentally  in  some  materials. 
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table  A1 

COMPUTER  CODE  FOR  CALCULATING  HUGONIOTS  FOR  POROUS  SOLIDS 

100  PROGRAM  HUG220  (  I NPUT fOUTPUT ) 

102  40  PRINT  600 
104  PRINT  650 

106C  TYPICAL  VALUES  ARE  0.239»  2.05»  0.737,  1.614,  0.62 

108  READ  660»C0*S.G0*A0*VS 

no  PRINT  660*  C0*S*G0*A0.VS 

112  VF  =  AO  *  VS 

114  PRINT  600 

116  F  =  40.0 

118  V2  =  VS  *  AO  /  (1.0  +  2.0  /GO) 

120  V3  =  (VS  -  V2)  /  F 
122  V  =  VS 

124  PRINT  500.VS*V2»V3 
126  PRINT  600 
128  PRINT  800 


130 

DO 

30  I  = 

1*20 

132 

PH 

=r 

CO  * 

CO  *  (VS  -  V)  / 

(VS  - 

S 

*  (VS  - 

V)  )**2 

134 

PP 

= 

PH  * 

(  1.0  -  GO  *  (VS 

-V)  ♦ 

0. 

5  /  V) 

136 

PP 

PP  / 

(1.0  -  0.5  *  GO 

*  (AO 

* 

VS  -  V) 

/  V) 

138 

U8 

= 

VF  * 

SQRT(PP  /  (VF  - 

V)  ) 

140 

U9 

= 

SORT 

(PP  *  ( VF  -  V) ) 

142 

El 

r 

1.195E4  ♦  U9**2 

144  PRINT  700.V.PH*PP*U8.U9*EI 

146  IF  (PP  .GE,  0.1)  GO  TO  40 

148  V  =  V  -  V3 

150  30  CONTINUE 

152  500  FORMAT  ( 3F10.4) 

154  600  FORMAT  ( IH  ) 

156  650  FORMAT! 1H.*TYPE  VALUES  OF  CO*  $*  GO*  AO  AND  VS*  FORMAT  F8.4») 
158  660  FORMAT  (5F8.4) 

160  700  FORMAT  ( 3F 1 0 . 4 * 2 F 1 0 • 3 , F 10 . 1  ) 

162  800  FORMAT  ( 5X , # V* * 9X , *PH« , 8X ,*PP* ♦ 8X * *08* * 8X * *09* * 8X * *E I*  ) 

164  END 
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table  A2 

OUTPUT  OF  CODE  LISTED  IN  TABLE  A1 
(TNT  AT  1.0  G/CC*  AO  =  1.614) 


TYPE  VALUES  C0»  S»  GO*  AO*  AND  VS 
0*239  2.050  0*737  1.614  0.62 

.2390  20500  .7370  1.6140  .6200 


.6200 

.2595 

.0088 

V 

PH 

PP 

U8 

U9 

El 


6200 

0.0000 

0.0000 

6112 

.0014 

.0018 

6025 

.0029 

.0038 

5937 

.0047 

.0062 

5849 

.0067 

.0088 

5762 

.0089 

.0119 

5674 

.0114 

.015  4 

5587 

.0143 

.0194 

5499 

.0177 

0241 

,5411 

.0214 

.0295 

,5324 

.0258 

.0359 

,5236 

.0309 

.0433 

,5148 

.0367 

*0521 

,5061 

.0436 

.0625 

,4973 

.0516 

.0748 

,4885 

.0611 

.0897 

,4798 

.0724 

.1077 

,4710 

.0860 

.1207 

,4623 

.1024 

.1569 

,4535 

.1225 

.1907 

TYPE  VALUES  OF  CO*  S*  GO*  AO* 


0.000 

0*000 

0*0 

.068 

.026 

8*4 

.098 

.039 

18*3 

.123 

.050 

30*0 

.146 

.061 

43*9 

.167 

.071 

60*3 

.189 

.082 

79*7 

.210 

.093 

102*6 

.231 

.104 

129*9 

.254 

.117 

162*2 

.277 

.130 

200*8 

.302 

.144 

246.9 

.328 

.159 

302.3 

.356 

.176 

369.1 

.386 

.194 

450*1 

*419 

.214 

549*0 

.455 

.237 

670*5 

.495 

.262 

821*1 

.540 

.291 

1009*3 

*591 

.323 

1247.0 

AND  VS 
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TABLE  A3 

COMPUTED  HUGONIOTS  FOR  TNT  AT  7  DENSITIES 


CO  =  0.239*  S  =  2.05*  G0=  0.737*  AO  =  1.5084*  VS  =  0.62 

DENSITY  =  1.07  G/CC 


y 

PH 

PP 

U8 

U9 

El 

.6200 

O.OOCO 

0.0000 

0.000 

0.000 

0.0 

.6108 

.0015 

.0018 

.070 

.024 

7.0 

.6016 

.0031 

.0039 

.101 

-036 

15.4 

.5924 

.0050 

.0062 

.126 

.046 

25.4 

.5832 

.0071 

.0089 

.149 

.056 

37.5 

.5740 

.0095 

.0120 

.171 

.066 

51.9 

.5648 

.0123 

.0156 

.192 

.076 

69.1 

.5556 

.0155 

.0198 

.213 

.087 

89.7 

.5464 

.0191 

.0246 

.235 

.098 

114.4 

.5372 

.0233 

.0303 

.258 

.110 

144.1 

.5280 

.0283 

.0369 

.282 

.123 

179.8 

.5)88 

.0340 

.0448 

.307 

.137 

223.0 

.5095 

.0407 

.0541 

.334 

.152 

275.4 

.5003 

.0487 

.0653 

.362 

.169 

339.4 

.4911 

.0581 

.0783 

.394 

.187 

417.9 

.4819 

.0695 

.0951 

.428 

.208 

515.0 

.4727 

.0831 

.1151 

.467 

.231 

636.0 

.4635 

.0998 

.1398 

.509 

.257 

788.3 

.4543 

.1204 

.1709 

.557 

.287 

981.9 

.4451 

.1461 

.2102 

.613 

.321 

1231.3 

=  0.239, 

S  =  2.05* 

G0=  0.737 

*  AO  =  1.2912 

*  VS  = 

0.62 

DENSITY  =  1. 

25  G/CC 

V 

PH 

PP 

U8 

U9 

El 

.6200 

0.0000 

0.0000 

0.000 

0.000 

0.0 

.6099 

.0016 

.0018 

.078 

.019 

4.1 

.5998 

.0035 

.0039 

.111 

.028 

9.3 

.5897 

.0056 

.0063 

.138 

.036 

15.9 

.5796 

.0080 

.0091 

.162 

.045 

24.0 

.5694 

.0108 

.0123 

.185 

.053 

34.0 

.5593 

.0141 

.0161 

.207 

.062 

46.4 

.5492 

.0179 

.0205 

.229 

.072 

61.7 

.5391 

.0224 

.0258 

.251 

.082 

80.5 

.5290 

.0277 

.0320 

.275 

.093 

103.7 

.5189 

.0339 

.0393 

.299 

.105 

132.4 

.5083 

.0413 

.0482 

.325 

.119 

168.0 

.4987 

.0503 

.0589 

.354 

.133 

212.5 

.4886 

.0611 

.0720 

.385 

.150 

268.4 

.4784 

.0743 

.0881 

.419 

.168 

339.0 

.4683 

.0907 

.1081 

.457 

.190 

429.2 

.4582 

.1111 

.1334 

.500 

.214 

545.7 

.4481 

.1371 

.1657 

.549 

.242 

698<.0 

.4380 

.1705 

.2078 

.606 

.274 

900.2 
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0*239* 


.6200 
•  6092 
.5986 
.5877 
.5770 
.5662 
.5555 
.5447 
.5340 
.5232 
.5124 
.5017 
.4909 
.4802 
.4694 
.4587 
.4479 
.4371 


0.239, 


.6200 
.6090 
.5980 
.5871 
.  5761 
.5651 
.5541 
.  543  1 
.  5322 
.5212 
.5102 
.4992 
.4883 
.4773 
.4663 
.4553 
.4443 
.4334 


TABLE  A3  (CONTINUED) 


=  2.05,  GO 

=  0.737, 

AO  =  1.1366 

,  VS  = 

0.62 

DENS  I 

TY  =  1.42 

G/CC 

PH 

PP 

U8 

U9 

El 

0.0000 

0.0000 

0.000 

0.000 

0.0 

.0017 

.0018 

.097 

.013 

2.1 

.0037 

.0039 

.135 

.020 

5.0 

.0060 

.0064 

.164 

.027 

8.9 

.0087 

.0092 

.189 

.034 

14.0 

.0118 

.0125 

.212 

.042 

20.7 

.0155 

.0165 

.234 

.050 

29.4 

.0198 

.0211 

.256 

.058 

40.4 

.0250 

.0266 

.278 

.067 

54.3 

.0311 

.0332 

.302 

.078 

72.1 

.0385 

.0412 

.326 

.089 

94.7 

.0474 

.0509 

.353 

.102 

123.5 

.0584 

.0628 

.382 

.116 

160.4 

.0719 

.0775 

.414 

.132 

208.0 

.0888 

.0960 

.450 

.150 

269.9 

.1102 

.1195 

.491 

.171 

351.3 

.1377 

.1499 

.538 

.196 

459.8 

.1738 

.1898 

.594 

.225 

606.9 

=  2.05, 

G0=  0.''37 

,  AO  =  1.0832 

,  VS  ^ 

0.62 

DENSITY  =  1. 

49  G/CC 

PH 

PP 

U8 

U9 

El 

o.onon 

0.0000 

0.000 

0.000 

0.0 

.0018 

.0018 

.114 

.011 

1.4 

.0038 

.0039 

.155 

.017 

3.4 

.0062 

.0064 

.184 

.023 

6.4 

.0089 

.0092 

.209 

.030 

10.6 

.0122 

.0126 

.231 

.037 

16.1 

.0160 

.0166 

.252 

.044 

23.3 

.0205 

.0213 

.274 

.052 

32.7 

.0259 

.0260 

.295 

.061 

44.9 

.0324 

.0337 

.318 

.071 

60.6 

.0402 

.0419 

.342 

.082 

80.8 

.0497 

.0519 

.369 

.095 

106.9 

.0615 

.0642 

.398 

.109 

140.7 

.0761 

.0796 

.430 

.124 

184.9 

.0944 

.0990 

.466 

.143 

243.0 

.1180 

.1239 

.508 

.164 

320.3 

.1486 

.1564 

.557 

.189 

424.7 

.1892 

.1996 

.615 

.218 

568.2 
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TABLE  A3  ^CONTINUED) 


0.239f 

S  =  2.05. 

G0=  0.737 

.  AO  =  1.0215 

*  VS  = 

0.62 

DENSITY  =  1. 

58  G/CC 

PH 

PP 

U8 

U9 

El 

.6200 

0.0000 

0.0000 

0.000 

0.000 

C.O 

•  6068 

.0018 

.0018 

.172 

.007 

.5 

.5975 

.0039 

.0039 

.210 

.012 

1.7 

.5863 

.0063 

.0064 

.234 

.017 

3.6 

.5751 

.0092 

.0093 

.253 

.023 

6.5 

.5638 

.0126 

.0127 

.271 

.030 

10.6 

.5526 

.0166 

.0167 

•  288 

.037 

16.2 

.5413 

.0213 

.0216 

.307 

.045 

23.7 

.5301 

.0270 

.0273 

.326 

.053 

33.7 

.5189 

.0339 

.0343 

.347 

.063 

46.9 

.5076 

.0423 

.0427 

.369 

.073 

64.2 

.4964 

.0525 

.0531 

.394 

.085 

86.9 

.4852 

.0652 

.0660 

.423 

.099 

116.8 

.4739 

.0812 

.0822 

.455 

.114 

156.5 

.4627 

.1015 

.1028 

.491 

.132 

209.6 

.4515 

.1278 

.1294 

.534 

.153 

281.3 

.4402 

.1624 

.1646 

.585 

.178 

379.8 

.4290 

.2092 

.2121 

.645 

.208 

517.8 

0.239. 

S  =  2.05. 

G0=  0.737 

.  AO  =  1.0088 

*  VS  = 

0.62 

density  =  1. 

60  G/CC 

PH 

PP 

U8 

U9 

El 

.6200 

0.0000 

0.0000 

0.000 

0.000 

0.0 

.6087 

.0018 

.0018 

.206 

.006 

.4 

.5974 

.0039 

.0039 

.234 

.010 

1.3 

.5861 

.0064 

.0064 

.252 

.016 

3.0 

.5748 

.0093 

.0093 

.268 

.022 

5.6 

.5636 

.0127 

.0127 

.264 

.028 

9.4 

.5523 

.0167 

.0168 

.299 

.035 

14.7 

.5410 

.0215 

.0216 

.316 

.043 

21.8 

.5297 

.0273 

.0274 

.335 

.051 

31.4 

.5184 

.0342 

.0344 

.354 

.061 

44.0 

.5071 

.0427 

.0429 

.377 

.071 

60.7 

.4958 

.0531 

.0534 

.401 

.083 

82.7 

.4845 

.0661 

.0664 

.429 

.097 

111.8 

.4732 

.0823 

.0827 

•  461 

.112 

150.5 

.4619 

.1030 

.1036 

.498 

.130 

202.4 

.4507 

.1299 

.1306 

.541 

.151 

272.8 

.4394 

.1655 

.16^4 

.591 

.176 

370.0 

.4281 

.2136 

.2148 

.652 

.206 

506.6 
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TABLE  A3  (CONTINUED) 

0.239f  S  =  2.05*  G0=  0.737*  AO  =  1.0000*  VS  =  0.62 


DENSITY  ^  1. 

614  6/CC 

PH 

PP 

U8 

U9 

El 

.6200 

o.onoo 

0*0000 

0.000 

0.000 

0.0 

.6087 

.0018 

.0018 

.248 

.005 

.2 

.5973 

.0039 

.0039 

.258 

.009 

1.1 

.5860 

.0064 

.0064 

.269 

.015 

2.6 

.5747 

.0093 

.0093 

.281 

.021 

5.0 

.5634 

.0127 

.0127 

.294 

.027 

8.6 

.5520 

.0168 

.0168 

.308 

.034 

13.6 

.5407 

.0216 

.0216 

.324 

.041 

20.5 

.5294 

.0274 

.0274 

.341 

.050 

29.7 

.5181 

.0345 

.0345 

.361 

.059 

42.0 

.5067 

.0430 

.0430 

.382 

.070 

58.2 

,4954 

.0535 

.0535 

.406 

.082 

79.7 

.4841 

.0666 

.0666 

.434 

.095 

108.2 

.47:>'‘ 

.0831 

.0831 

.466 

.111 

146.2 

.  4o  14 

.1041 

.10^1 

.502 

.128 

197.3 

.4501 

.1314 

.1314 

.545 

.149 

266.9 

.4388 

.1676 

.1676 

.596 

.174 

363.0 

.4275 

.2167 

.2167 

.658 

.204 

498.6 
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TABLE  A4 

COMPUTED  HUGONIOTS  FOR  5  EXPLOSIVES  USING  RUSSIAN  DATA 


EXPLOSIVE  -  NB  AT  1.00  G/CC 


=  0.170* 

S  =  1.85* 

GO  =1.50* 

AO  =  1.63* 

VS  = 

0.6135 

V 

PH 

PP 

U8 

U9 

El 

.6135 

0.0000 

0.0000 

0.000 

0.000 

0.0 

.6043 

.0008 

.0015 

.061 

.024 

6.9 

.5950 

.0016 

.0032 

.089 

.036 

15.4 

.5858 

.0025 

.0052 

.112 

.046 

25.8 

.5765 

.0036 

.0076 

.134 

.057 

38.6 

.5673 

.0048 

.0105 

.156 

.067 

54.4 

.5580 

.0061 

.0140 

.178 

.079 

74.0 

.5488 

.0077 

.0182 

.201 

.091 

98.4 

.5395 

.0094 

.0235 

.226 

.104 

129.1 

.5303 

.0114 

.0299 

.252 

.119 

168.0 

.5210 

.0137 

.0381 

.282 

.135 

218.1 

.5118 

.0162 

.0486 

.315 

.154 

283.5 

.5025 

.0192 

.0623 

.354 

.176 

370.5 

.4933 

.0227 

.0808 

.399 

.202 

489.4 

.4841 

.0267 

.1066 

.455 

.235 

657.2 

.4748 

.0315 

.1442 

.524 

.275 

905.0 

.4656 

.0370 

.2029 

.616 

.329 

1295.6 

EXPLOSIVE  - 

RDX  AT  l.OC 

G/CC 

=  0.287* 

S  =  1.61* 

GO  =2.€0* 

o 

00 

. 

t— 1 

II 

o 

< 

VS  = 

0.5556 

V 

PH 

PP 

U8 

U9 

El 

.5556 

0.0000 

0.0000 

0.000 

0.000 

0.0 

.5561 

-.0001 

.0034 

.087 

.039 

18.0 

.5566 

-.0003 

.0072 

.127 

.056 

37.9 

.5570 

-.0004 

.0113 

.160 

.071 

60.1 

.5575 

-.0005 

.0161 

.191 

.084 

84.9 

.5580 

-.0006 

.0214 

.220 

.097 

113.0 

.5585 

-.0008 

.0275 

.250 

.110 

145.1 

.5590 

-.0009 

.0345 

.280 

.123 

181.9 

.5595 

-.0010 

.0427 

.311 

.137 

224.8 

.5599 

-.0011 

.0524 

.345 

.152 

275.4 

.5604 

-.0013 

.0639 

.381 

.168 

335.9 

•  5609 

-.0014 

.0781 

.422 

.185 

409.8 

.5614 

-.0015 

.0958 

.467 

.205 

502.0 

.5619 

-.0016 

.1185 

.520 

.228 

620.4 

.5624 

-.0017 

.1487 

.583 

.255 

777.9 

.5628 

-.0019 

.1911 

.661 

.289 

998.3 
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CO  = 


CO  = 


TABLE  A4  (CONTINUED) 


EXPLOSIVE  -  TETRYL  AT  0.86  G/CC 


.217, 

S  =  1.91, 

GO  =1.65, 

AO  =  2.012 

,  VS  = 

0.5700 

PH 

PP 

U8 

U9 

El 

5700 

0.0000 

0.0000 

0.000 

0.000 

0.0 

5754 

.0004 

.0024 

.074 

.037 

16.6 

5728 

.0008 

.0051 

.108 

.055 

35.6 

5702 

.0012 

.0081 

.136 

.069 

57.4 

5675 

.0016 

.0116 

.162 

.083 

82.4 

5649 

.0020 

.0156 

.188 

.097 

111.4 

5623 

.0025 

.0202 

.213 

.110 

145.2 

5597 

.0029 

.0257 

.240 

.124 

105.0 

5571 

.0034 

.0321 

•  268 

.139 

232.2 

5545 

.0039 

.0390 

.297 

.156 

209.0 

5519 

.0044 

.0491 

.330 

.173 

358.3 

5492 

.0049 

.0606 

.365 

.193 

444.2 

5466 

.0055 

.0751 

.406 

.215 

553.0 

5440 

.0061 

.0939 

.453 

.241 

694.6 

5414 

.0067 

.1192 

.509 

.272 

885.5 

5308 

.0073 

.1549 

.579 

.311 

1155.2 

5362 

.0079 

.2087 

.671 

.362 

1563.0 

EXPLOSIVE  - 

PETN  AT  0. 

82  G/CC 

.242, 

S  =  1.91, 

GO  =1.70» 

AO  =  2.159 

,  VS  = 

0.5650 

PH 

PP 

U8 

U9 

El 

5650 

0.0000 

0.0000 

0.000 

0.000 

0.0 

5640 

.0000 

.0030 

.082 

.044 

23.1 

5645 

.0001 

.0062 

.119 

•  064 

48.8 

5643 

.0001 

.0099 

.150 

.081 

77.6 

5641 

.0002 

.0141 

.179 

.096 

110.1 

5639 

.0002 

.0188 

.206 

.111 

147.0 

5636 

.0003 

.0241 

.234 

.126 

189.2 

5634 

.0003 

.0303 

.262 

.141 

238.0 

5632 

.0003 

.0376 

.292 

.157 

295.0 

5630 

.0004 

.0462 

.323 

.174 

362.4 

5627 

.0004 

.0565 

,358 

.193 

4 '.3, 5 

5625 

.0005 

.06C'l 

.395 

.213 

542.6 

5623 

.0005 

.0848 

.438 

.236 

6S6.7 

5621 

.0006 

.1051 

.488 

.263 

8/6.4 

5618 

.0006 

.1322 

.547 

.295 

1039.5 

5616 

.0006 

.1701 

.620 

.335 

1338.0 

5614 

.0007 

.2270 

.716 

.387 

1786.1 

AlO 


UNCLASSIFIED 


UNCLASSIFIED 
NOLTR  74-213 


TABLE  A4  (CONTINUED) 


EXPLOSIVE  -  AN  AT  0.86  G/CC 


0.220* 

S  =  1.96* 

GO  =0.90* 

AO  =  2.012 

*  VS  = 

0.5780 

V 

PH 

PP 

U8 

U9 

El 

2200  1 

.9600  .9000  2.0120 

.5780 

.5780 

0.0000 

0.0000 

0.000 

0.000 

0.0 

.5671 

.0017 

.0032 

.085 

.044 

22.7 

.5563 

.0037 

.0071 

.126 

.065 

51.3 

.5454 

.0060 

.0116 

.161 

.085 

87.3 

.5346 

.0086 

.0177 

.195 

.105 

132.8 

.5237 

.0118 

.0250 

.230 

.126 

190.8 

.5129 

.0155 

.0341 

.266 

.149 

265.0 

.5020 

.0200 

.0457 

.306 

.174 

360.7 

.4912 

.0253 

.0605 

.349 

.202 

485.6 

.4803 

.0316 

.0798 

.398 

.233 

650.7 

.4695 

.0394 

.1052 

.453 

.270 

872.1 

.4586 

.0486 

.1396 

.518 

.314 

1175.0 

.4477 

.0605 

.1871 

.595 

.366 

1599.0 

All 
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APPENDIX  B 

ESTIMATED  ERRORS  IN  THE  PRESSURE  AND  THE  SPECIFIC  VOLUME 

As  shown  in  the  section  on  Inipedance  Matching ,  the  pressure  induced 
in  the  porous  TNT  is  obtained  graphically.  Analysis  of  the  error  in 
such  a  process  is  difficult.  A  mathematical  model  would  simplify  the 
estimation  of  these  errors.  The  impedance  matching  process  can  be 
described  mathematically  for  the  case  in  which  the  relief  path  in  the 
P,u  plane  for  the  PMMA  is  a  mirror  image  of  the  Hugoniot  curve.  Thus, 
to  estimate  errors  we  must  ignore  the  viscoelastic  effects  on  the 
unloading  path  for  the  PMMA.  Because  of  the  peculiar  U,u  relation  for 
PMMA  at  pressures  below  about  20  kbar,  the  following  analysis  does  not 
apply  for  low  pressure  work. 

In  the  following,  it  is  assumed  that  the  relation  between  shock 
velocity,  U  and  particle  velocity,  u,  is  linear,  so  that 

U  =  Co  +  Su,  (Bl) 

where  velocities  are  in  mm/ysec.  This  relation  holds  for  pressures 
greater  than  21.7  kbar  in  PMMA^.  Note  that  when  we  have  this  pressure 
in  the  PMMA,  we  will  have  a  considerably  lower  pressure  in  the 
approximately  1.0  g/cc  TNT.  Using  the  relation 

P  =  lOpo  u  U  (B2) 

where  po  is  the  density  (1.185  g/cc)  of  the  unshocked  PMMA,  and  P  is 
the  pressure  in  kbars,  we  have 

P  =  10  Po  u  (cq  +Su)  =Au+B  u^.  (B3) 

The  coefficients  in  equation  Bl  have  the  values  2.56  ±  0.01  mm/psec 
and  1.595  ±  0.01  respectively^,  where  the  uncertainty  is  a  standard 
deviation  as  given  by  a  least  squares  fitting  algorithm.  Thus,  the 

coefficients  in  equation  B3,  A  and  B,  are  30.3  ±  0.1  and  18.9  ±  0.11 

respectively. 

The  mirror  image  of  the  Hugoniot  curve  on  the  line  u  =  u^,  see 
Figure  5  in  the  main  text,  is  described  by 

P  =  A{2  u^  -u)  +  B(2  -u)^.  (B4) 

The  pressure  corresponding  to  the  particle  velocity,  Ui,  is  Pi.  As 

shown  in  the  section  on  Impedance  Matching,  we  must  locate  the  inter¬ 
section  of  a  curve  described  by  Equation  B4  with  a  line  u  =  Ug,  where 
Ug  is  the  quantity  measured  by  the  EMV  gage  at  the  interface  between 
the  PMMA  and  the  porous  TNT.  The  pressure  at  the  point  where  these 
curves  cross  is  given  by 

Pg  =  A(2  ui  -Ug)  +  B(2  Ui  -Ug)^.  (B5) 
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The  specific  volume  is  given  by 

Ve  =  “  10  “e^/Pe  (B6) 

where  Vof  is  the  initial  volume  of  the  porous  material.  We  are  con¬ 
cerned  with  the  effect  on  the  values  of  and  Vq  of  errors  in  the 
quantities  A,  B,  Ug.  We  will  assume  that  Vof  is  without  error. 

Differentiating  the  above  two  equations  gives 

dPg  =  dA(2ui-Ug)  +  dB  (2ui-Ug)^  +  dU]^  ( 2A+4B  ( 2ui-Ug)  ) 

-dUe  (A+2B(2ui-Ue) ) ,  (B7) 

and 

dVg  =  10Ue(Ue  dPe-2PedUe) /Pe^*  (B8) 

Note  that  dPg  depends  on  four  terms,  each  of  which  contains  a  differen¬ 
tial  of  a  variable  or  parameter  which  may  be  in  error.  As  is  conven¬ 
tional  in  error  analysis,  we  replace  the  differentials  by  finite 
differences,  and  add  the  absolute  values  of  the  individual  terms, 
giving, 

APg  =  I  AA  (2U]_-Ug)  I  +  I  AB  (2U]_-Ue)  ^  I  +  |  Au]_  (2A+4B  (2u]_-Ug) )  | 

+  I AUg (A+2B (2ui-Ug) ) I  (B9) 

and 

AVg  =  lOUedugAPel  +  |  2PeAUe  | ) /Pe^  .  (BIO) 

Equations  such  as  B8  and  B9  usually  overestimate  the  uncertainty  in 
the  error.  Another  way  of  combining  the  terms  is  to  take  the  square 
root  of  the  sum  of  the  squares  of  the  individual  terms.  In  the 
following,  such  estimates  will  be  indicated  by  placing  a  bar  over  the 
symbol  representing  the  error,  i.e.,  the  estimated  errors  in  Pg  and 
Vg  are  APg  and  Avg  when  this  method  is  used. 

In  order  to  evaluate  the  errors  as  given  by 'Equations  (B9)  and  (BIC) 
we  need  values  of  the  estimated  errors  in  A,  B,  uj ,  and  Ue.  The 
estimated  errors  in  A  and  B  are  given  above--the  values  are  from  the 
computer  code  that  was  used  to  fit  the  linear  expression,  Equation  Bl, 
to  the  Hugoniot  data  of  PMMA.  The  error  in  the  particle  velocity  as 
measured  by  the  EMV  gage  has  been  given  before^.  In  that  work,  it 
was  shown  that  the  error  in  the  particle  velocity  was  made  up  of  7 
components.  The  most  important  components  were  the  error  involved  in 
reading  the  calibration  marks  on  the  oscilloscope  record  and  in  reading 
the  oscilloscope  trace  itself.  The  total  error  was  taken  to  be  the 
sum  of  the  absolute  values  of  the  seven  terms.  Using  the  nomenclature 
adopted  above  we  are  presenting  Au.  The  results  are  given  in  column  2 
in  Table  Bl.  Note  that  there  is  some  dependency  of  the  error  on  the 
magnitude  of  the  measured  quantity.  The  variation  is  small,  however. 
The  reason  for  this  is  that  the  sensitivity  of  the  oscilloscope  is 
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TABLE  B1 

ESTIMATED  ERRORS  IN  PARTICLE  VELOCITY  AS  MEASURED  WITH  EMV  GAGE 


Particle  Velocity,  u 
mm/ysec 

Au  X  100% 
u 

Au  X  100% 
u 

0.1 

6.2 

3.1 

1.0 

5.7 

3.0 

2.0 

5.6 

3.0 

changed  so  that  approximately  the  same  deflection  is  obtained  for  all 
velocities.  In  Reference  3  it  was  noted  that  the  variation,  such  as 
it  is,  depends  entirely  on  the  3  digit  voltmeter  which  is  used  to 
calibrate  the  oscilloscope.  The  estimate  of  the  error  in  the  third 
column  was  obtained  by  taking  the  square  root  of  the  sum  of  the 
squares  of  the  seven  terms  contributing  to  the  error  in  the  particle 
velocity.  Again,  using  the  nomenclature  developed  above,  this  quantity 
is  identified  as  Au.  The  third  column  of  Table  B1  consists  of  values 
arrived  at  by  this  latter  process.  For  these  estimated  errors,  the 
variation  with  the  magnitude  of  the  particle  velocity  is  negligible. 
Also,  the  estimated  error  is  now  considerably  smaller,  consistent  with 
the  expectation  that  some  of  the  random  errors  should  compensate  for 
others.  lii  the  following  we  use  ±3%  as  the  error  in  the  values  of 
both  ui  and  u^  and  we  proceed  to  evaluate  Equations  (B9)  and  (BIO) . 

This  means  we  have  changed  Au]^  and  AUg  in  Equations  (B9)  and  (BIO)  to 
AU]^  and  AUg. 

We  have  chosen  four  pairs  of  values  of  u^  and  u^  from  the  raw 
data  of  Table  1  of  the  main  text.  These  are  listed  in  Table  B2,  along 
with  the  corresponding  values  of  and  Pg .  Note  that  P]_  varies  from 
20.7  to  57.1  kbars,  while  Pg  varies  from  12,5  to  43.0  kbars,  which 
covers  the  upper  range  of  our  data.  The  low  pressure  point,  for  which 
P]^  is  20.7  kbar,  requires  a  small  extrapolation  of  the  linear  U,u 
relation.  Equation  (Bl) ;  this  introduces  a  negligibly  small  error  in 
the  value  of  P]^.  The  values,  in  kbars,  of  the  individual  terms  of 
Equation  (B9)  are  given  in  the  columns  headed  Ei  through  E^.  Note 
that  both  E^  and  E2  are  negligibly  small.  In  the  last  column  we  show 
the  square  root  of  the  sum  of  the  squares  of  the  4  components  of  the 
errors,  i.e.,  APg.  These  errors,  which  are  in  kbars,  must  be  compared 
with  the  values  of  Pg  which  are  given  in  column  5  of  Table  B2.  We 
note  that  the  magnitude  of  the  error  varies  with  the  value  of  P^.  But 
the  error  on  a  percentage  basis,  see  the  last  column  of  the  table,  does 
not  vary  significantly  with  the  pressure.  Apparently  our  error  in  the 
pressure  in  the  test  explosive  is  12  to  13%  over  a  range  of  pressure 
from  12.5  to  43.0  kbars  in  the  low  aensity  TNT. 

The  error  in  the  specific  volume,  Vg,  of  the  shock  compressed 
sample  is  computed  from  Equation  (BIO) in  which  APg  and  Ue  are  replaced 
^respectively  by  APg  and  Aug.  We  have  used  the  same  pairs  of  values  of 
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and  Ug  that  we  used  for  computing  the  error  in  the  pressure,  see 
Table  B2.  Equation  (B6)  is  used  to  obtain  the  values  of  Vg,  see  the 
first  column  of  Table  B3.  Values  of  Pg  and  APe  are  from  Table  B2.  In 
this  case,  the  error  has  two  components ,  labeled  Fi  and  F2  in  Table  B3. 
The  first  component,  that  involving  AP^,  dominates  the  second. 

Combining  the  errors  by  squaring,  adding,  and  then  taking  the  square 
root  gives  the  expected  error  in  volume,  Ave,  see  the  next  to  last 
column  in  Table  B3.  In  the  last  column,  the  errors  are  converted  to 
percent  relative  errors.  We  have  estimated  errors  ranging  from  about 
7%  to  13%  in  the  data  of  Table  1,  i.e.,  the  raw  data. 
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APPENDIX  C 

MELTING  OF  TNT  UNDER  SHOCK  LOADING 

Shock  induced  phase  changes,  such  as  melting,  in  materials  are 
common.  Boyle,  Smothers,  and  Ervin^^  indicate  that  melting  may  occur 
in  cast  TNT,  =  1.62  g/cc,  at  around  32  kbar.  In  the  section 
entitled  "Comparison  of  Calculated  and  Experimental  Hugoniot  for  Porous 
Explosives,  Data  from  NOL  and  Other  Sources",  the  possibility  of 
melting  of  0.98  g/cc  TNT  was  raised  as  one  explanation  of  the  data 
points  at  P  -  8.5  kbar,  v  =  0.795  cc/g,  see  Figure  Cl.  In  this 
appendix,  this  possibility  will  be  discussed  more  fully. 

If  the  TNT  is  melting  at  P  -  8.5  kbar,  v  =  0.795  cc/g  (or  close 
to  this  point) ,  the  porous  Hugoniot  would  have  to  be  based  on  the 
parameters  of  liquid  TNT.  These  are^^ 

Po  =  1.472  g/cc  (Cl) 

U  =  2.14  +  1.57  u  ,  T  =  81®C  (C2) 

where  U  and  u  are  in  mm/ysec.  Fq  is  assumed  to  be  the  same  as  that 
for  the  solid,  0.737.  The  resulting  Hugoniot  is  shown  in  Figure  Cl. 

As  is  evident  from  the  figure,  our  data  agrees  with  either  Hugoniot 
equally  well. 

In  shock  loading  studies,  the  usual  procedures  employed  to  detect 
a  phase  change  are  (1)  determine  if  there  is  a  break  in  the  Hugoniot, 
and  (2)  determine  if  the  initial  wave  splits  into  two  waves.  (An 
excellent  discussion  of  shock  induced  phase  changes  is  given  by  McQueen, 
Marsh,  Taylor,  Fritz,  and  Carter^**  .)  For  0.98  g/cc  TNT,  a  break  in 
the  Hugoniot  may  exist  at  the  above  two  mentioned  data  points. 

However,  these  data  may  also  be  due  to  crushup  (or  both) .  The 
temperature  rise  at  these  two  data  points  can  be  calculated  from  the 
following  equation 


AT  =  [1/2(P  +  Pq) (Vof  -  V)  -  H.F.l/c^  (C3) 

where  H.F.  is  the  latent  heat  of  fusion  and  is  equal  to  22.34  cal/g 
for  TNT  at  79®C^^ ,  and  c^  is  the  specific  heat  at  constant  volume  and 
is  equal  to  0.29  cal/g/®C  (see  calculation  of  values  of  Tq  in  Ref.  21). 

AT  for  these  two  data  points,  using  Equation  (C3) ,  are  60.9°C  and 
63.8®C.  The  initial  temperature  of  the  TNT  was  25°C  so  that  the  final 
temperatures  in  the  TNT  were  85.9®C  and  88.8"C.  The  melting  point  of 
TNT  at  ambient  pressure  is  81®C^®  Thus,  the  melting  of  the  TNT  is  a 
possibility  if  it  is  not  strongly  dependent  on  pressure. 

In  an  effort  to  resolve  this  question,  an  attempt  was  made  to 
detect  multiple  waves.  Several  experiments  were  performed  with  EMV 
gages  located  in  the  TNT  at  different  distances  from  the  PMMA/TNT 
interface.  The  length  of  the  P^‘MA  attenuator  was  50.8  mm.  The  inter¬ 
face  values  of  u^  and  P^  were  -  0.70  mm/ysec  and  -  12  kbar,  respectively. 
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This  is  well  above  the  values  of  the  data  which  may  indicate  a  phase 
change  (i.e.,  at  -  8.5  kbar)  but  below  the  Rayleigh  line  drawn  through 
these  data  at  -  8.5  kbar  from  Vof .  Thus,  if  a  phase  change  is 
present  in  the  vicinity  of  P  -  8.5  kbar,  u  0.44  mm/|isec,  and 

V  =  0.795  cc/g,  the  wave  should  split  for  these  initial  interface 
values. 

The  u,t  profiles  generated  by  the  EMV  gage  (with  mica  insulated 
leads)  located  in  the  TNT  are  shown  in  Figure  C2.  The  number  on  each 
curve  indicates  the  thickness  in  mm  of  the  TNT  between  the  PMMA/TNT 
interface  and  and  the  gage.  One  u,t  profile  from  a  gage  located  at 
9  mm  had  a  spurious  electrical  signal  superimposed  on  it.  This  profile 
was  smoothed  by  drawing  a  curve  through  the  noise.  This  is  the  only 
profile  which  showed  a  hint  of  wave  splitting.  The  other  9  mm  profile 
(and  an  additional  9  mm  profile  not  shown)  have  a  faster  risetime  and 
reach  a  higher  particle  velocity  than  the  interface  (0  mm)  result. 
These  two  facts  indicate  that  buildup  to  detonation  is  occurring  and 
eliminates  any  possibility  of  seeing  wave  splitting.  Experiments  at 
X  >  9  mm  also  show  buildup.  Similar  experiments  were  performed  in 
1.30  g/cc  TNT  but  buildup  to  detonation  was  also  present,  see 
Appendix  E.  Because  buildup  to  detonation  occurs  so  rapidly  in  these 
experiments  we  have  not  been  able  to  confirm  or  reject  melting  of  the 
TNT  at  ~  8.5  kbar. 
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FIG.  C2  u,t  PROFILES  AT  VARIOUS  LOCATIONS  IN  0,98  g/cc  TNT 
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APPENDIX  D 

SOME  UNREACTED  HUGONIOT  DATA  FOR  TNT  AT  1.30  G/CC 

In  the  investigation  of  melting  of  TNT  (see  Appendix  C) ,  some 
experiments  were  performed  in  1.30  g/cc  TNT.  These  experiments  had 
one  EMV  gage  located  at  the  PMMA/TNT  interface  and  one  EMV  gage  located 
in  the  TNT.  The  gages  at  the  interface  yielded  unieacted  Hugoniot 
data. 


The  u,t  profiles  obtained  in  1.30  g/cc  TNT  showed  an  overshoot  in 
u  in  the  first  150  ns,  see  Figure  Dl.  This  overshoot  was  seen  in 
PMMA^  and  is  thought  to  be  due  to  a  piezoelectric  effect.  As  was  done 
in  PMMA^,  the  intercept  particle  velocity,  u^,  was  obtained  by  ignoring 
this  overshoot  and  beginning  the  extrapolation  at  ~  150  ns.  The 
resulting  u-  and  P^  (obtained  by  impedance  matching)  are  listed  in 
Table  Dl  and  plotted  in  Figure  D2.  Note  that  the  scale  of  Figure  D2 
is  -  twice  as  large  as  earlier  P,u  plots.  Because  these  data  were 
close  to  the  PMMA  Hugoniot  for  forward  facing  shock  (curve  1  of 
Figure  5),  the  correction  due  to  viscoelasticity  is  small  (-  0.2  kbar) 
compared  to  the  estimated  error  in  pressure  (2-3  kbar)  and  has  been 
ignored.  Also  plotted  in  Figure  D2  is  the  computed  Hugoniot  for 
1.30  g/cc  TNT  obtained  from  Herrmann's  model The  agreement  between 
the  experimental  data  and  the  calculated  Hugoniot  is  good. 
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TABLE  D1 


HUGONIOT  DATA  FOR  1.30  G/CC  TNT 


U 

V 

Pi 

X 

mm/ysec 

kbar 

mm/ysec 

cc/g 

Shot  No. 

iran/psec 

kbar 

mm 

0.581 

16.6 

2.198 

.566 

344 

.507 

20.3 

50.8 

0.609 

16.6 

2.097 

.546 

343 

.517 

20.7 

50.8 

0.730 

27.0 

2.845 

.571 

342 

.682 

29.4 

43.4 

D2 

UNCLASSIFIED 


UNCLASSIFIED 
NOLTR  74-213 


APPENDIX  E 

BUILDUP  OF  PARTICLE  VELOCITY  IN  LOW  DENSITY  TNT 

When  an  explosive  is  subjected  to  a  shock  stimulus  whose  impulse 
is  sufficiently  large  to  cause  eventual  detonation  but  not  great 
enough  to  cause  immediate  detonation,  the  particle  velocity,  shock 
velocity,  and  pressure  will  build  up  from  their  initial  values  to  their 
steady  state  detonation  values.  The  buildup  of  particle  velocity  in 
low  density  TNT  as  a  function  of  Z,  the  distance  from  the  inert/TNT 
interface  to  the  gage,  was  obtained  by  using  the  experimental  config¬ 
uration  shown  in  Figure  1  but  with  the  EMV  gages  placed  in  the  INT  at 
various  locations.  The  buildup  of  particle  velocity  was  investigated 
at  two  different  initial  pressures,  P^,  in  the  TNT,  4C  kbar  and 
-  12  kbar.  The  initial  density  was  0.98  g/cc  ±  0.03  g/cc*.  The 
quantity  of  interest  is  the  initial  value  of  the  particle  velocity,  u^. 
That  is,  the  velocity  at  zero  time  of  the  gage  record;  this  is  not  to 
be  confused  with  Uqj  which  occurs  at  some  time  t  >  o.  This  velocity 
is  obtained  by  fitting  the  u,t  result  to  a  straight  line  by  the  method 
of  least  squares.  The  results  are  presented  in  Table  El.  Some  of  the 
u,t  curves  obtained  at  various  values  of  i  for  the  same  initial  pressure 
in  the  TNT,  -  12  kbar,  are  plotted  in  Figure  El.  Zero  time  on  this 
figure  is  the  zero  time  of  the  gage  record.  Note  that  as  i  increases 
the  peak  of  the  particle  velocity  profile  increases  while  the  slope 
of  the  u,t  curve  after  the  peak,  in  most  cases,  becomes  more  negative. 
Reaction  will  cause  this  change  in  slope  as  will  the  gage  leads  being 
shorted  out  (see 'Effects  of  Reaction  on  Gage  Output").  The  EMV  gages 
used  in  the  shots  presented  in  Figure  El  did  not  have  their  leads 
encased  in  mica.  Thus,  the  rapid  falloff  of  the  u,t  curves  is  probably 
due  to  both  effects. 

Dremin,  Koldunov,  and  Shvedov^®  have  investigated  the  buildup  of 
both  particle  velocity  and  shock  velocity  as  a  function  of  I  in 
1.0  g/cc  and  0.78  g/cc  TNT  and  0.82  g/cc  tetryl  using  the  EMV  gage. 

At  1.0  g/cc  Dremin  investigated  two  different  particle  sizes  of  TNT, 

0.1  mm  and  0.37  -  1.0  mm.  Since  the  large  particle  size  used  by 
Dremin  is  closer  to  the  particle  size  of  the  TNT  used  here,  only  those 
data  will  be  compared  with  ours.  Instead  of  PMMA  as  the  attenuator, 
Dremin  used  paraffin.  The  diameter  of  their  charges  was  60  mm  while 
the  diameter  of  the  charges  used  in  this  study  was  50.8  mm.  The 
initial  pressures  in  the  Russian  work  were  8,  11.5,  and  25  kbar. 

Dremin  does  not  give  any  details  on  the  type  or  density  of  high  explo¬ 
sive  booster  he  used.  A  different  booster  system  will  cause  different 
buildup  results.  Results  for  the  two  higher  pressures  (open  symbols) 
are  compared  with  our  results  (solid  symbols)  in  Figure  E2.  The 
Russian  results  were  obtained  from  Figure  2  of  Reference  (20) .  At 
~  12  kbar,  our  results  and  Dremin 's  agree  well.  The  buildup  in 
particle  velocity  when  Pg  =  25  kbar  obtained  by  Dremin  appears  to  be 
more  rapid  than  the  buildup  obtained  in  this  study  when 

*Some  data  on  the  buildup  of  particle  velocity  in  1.30  g/cc  TNT  were 
also  obtained.  These  data  will  be  found  at  the  end  of  Table  El. 
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Pg  40  kbar.  This  may  be  due  to  the  difference  in  charge  diameter, 
particle  size,  and  booster.  The  largest  values  of  Ue  obtained  by 
Dremin  is  2.3  mm/psec  while  the  largest  values  observed  here  is 

2.2  mm/psec.  Again,  this  is  probably  due  to  the  difference  in  charge 
diameter,  particle  size,  and  booster. 

Dremin  also  obtained  the  relation  between  the  shock  velocity,  U, 
and  Ug  in  the  partially  reacting  TNT  at  1.0  g/cc.  His  result  is 

U  =  (0.30  +  1.85  u  )  mm/psec.  (El) 

Equation  (El)  is  practically  identical  with  our  Equations  (30)  and  (31) 
for  non- reacting  TNT  at  0.98  g/cc. 
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TABLE  El 


PARTICLE 

VELOCITY 

BUILDUP  DATA  FOR 

TNT 

£ 

ue 

Po 

Length  of  PMMA 

mm 

mm/psec 

g/cc 

Shot  No. 

X,  mm 

Pe  ^ 

40  kbar 

0.0 

1.34 

1.0 

291 

25.4 

6.5 

1.57 

1.0 

325A 

25.4 

25.0 

2.10 

1.0 

321 

25.4 

50.0 

2.19 

1.0 

322 

25.4 

50.0 

2.18 

1.0 

323 

25.4 

77.0 

- * 

1.0 

324 

25.4 

77.0 

- * 

1.0 

329 

25.4 

Pe  ' 

12  kbar 

0.0 

0.68 

1.02 

251 

50.8 

0.0 

0.73 

0.97 

253 

50.8 

3.2 

0.80 

0.98 

336 

46.8 

6.4 

0.90 

0.96 

279 

50.8 

9.0 

0.86 

0.98 

338™ 

50.8 

9.0 

0.80 

0.98 

339™ 

50.8 

12.0 

mm  — 

0.98 

338™ 

50.8 

12.4 

1.00 

0.96 

280 

50.8 

20.4 

1.21 

0.99 

281B 

50.8 

20.4 

1.27 

0.96 

;:85B 

50.8 

Pe  ' 

27  kbar 

0.0 

.730 

1.30 

342 

43.4 

6.4 

1.00 

1.30 

342 

43.4 

Pe  - 

16.6 

kbar 

0.0 

.58 

1.30 

344 

50.8 

0.0 

.61 

1.30 

343 

50.8 

6.4 

-.67 

1.30 

344 

50.8 

12.4 

-.80 

1.30 

343 

50.8 

Noise  on 

record. 

Risetime 

was  same 

as  total  sweep  of 

record. 

Mica  covered  gage 

leads . 
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FIG.  El  COMPOSITE  PLOT  OF  u,t  CURVES  (50.8  mm  OF  PMMA) 


FIG.  E2  COMPARISON  OF  BUILDUP  DATA 


